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ABSTRACT 


For a long time now, the analysis of a three phase power 
system network is being done on the assumption of a balanced 
system. However, due to certain inherent unbalances present 
in the system, it is necessary to evaluate the performance 
of the system in phase frame of reference. Moreover, the 
continuous increase in demand for power has exposed the 
limitations of three phase power systems. This has led to 
the feasibility studies of multiphase systems. In the present 
work, the detailed modelling of different elements present 
in three phase and multiphase (six phase and twelve phase) 
systems has been done. The modelling has been employed to 
carry out the balanced as well as unbalanced load flow analysi 
of three phase and six phase systems. 
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CHAPTER 1 


INTRODUCTION 

With the ever increasing demand for electric power and 
consequently, with the growing size and complexities of the 
modern power systems, it has become necessary to search for 
some viable alternatives to meet such a demand which requires 
transfer of more and more power from remote power plants to 
far away load centres. One of the alternatives, being utilized 
at present, is the high voltage transmission. However, trans- 
mission at such high voltages is not free of demerits. It 
poses some problems, for example, strong electric field at 
the ground surface with possible biological effects, high 
audible and radio noise, visual pollution, increasingly 
difficult problems in acquiring new rights of way, etc. 

Although work is being done to tackle such problems through 
improved design methods, search for other alternatives to 
transfer more power is also going on. 

One of the most promising alternatives to meet the 
increasing demand of electric power is the use of high phase 
order or multiphase systems) especially multiphase trans- 
mission. The idea of multiphase systems was first introduced 
by Barnes and Barthold in 1972 and ever since this topic has 
been an area of active research for power system engineers 
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and scientists* Feasibility studies of multiphase systems 
have shown that they offer considerable advantages over their 
three phase counterparts in several aspects. 

It is essential to develop an accurate model of multi- 
phase systems in order to study their feasibility. This is 
achieved by representing the different components of multi- 
phase systems in phasor coordinates with the help of mathe~ 
matical equations. The representation thus obtained is known 
as mathematical modelling of the multiphase systems which is 
a simple and inexpensive means to analyse the systems. The 
mathematical model, thus developed, is used to carry out the 
analysis of multiphase systems such as load flow analysis, 
fault analysis and transient stability analysis, etc* 

1*1 HISTORICAL SURVEY 

The concept of multiphase transmission was first 
introduced by Barnes and Barthold in 1972 [l]. This led J.R. 
Stewart ©t al [2,3] to carry out the feasibility studies of 
multiphase systems based on steady state considerations and 
their overvoltages and insulation requirements. Venkata et al 
, [5,6] pursued the topic by conducting preliminary investigation 
of power transfer capability, load flow, stability and relia- 
bility aspects of six phase systems. Grant et al [4] have 

K Confirmed analytical predictions of electrical and mechanical 

' ' 

«|$*|s^haviour of multiphase systems based on the studies of six 
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and twelve phase test systems, Laughton et al [7,12,15,16] 
have developed mathematical models of three phase transformers 
and machines in the phase frame of reference and have con- 
ducted the load flow and fault analysis of three phase systems, 
Roy et al [17] have also carried out the load flow analysis 
of three phase systems. The usual short circuit analysis of 
six phase transmission system was presented by Bhatt et al [13] 
and Nanda et al [19], Venkata et al [21] have discussed the 
various types of faults and their significant combinations 
likely to occur on a six phase line. The analytical expre- 
ssions for all the 23 significant faults have been developed. 
The effect of converting a double circuit three-phase line to 
a six-phase line was studied by L.P. Singh et al [10], Furhter, 
theoretical work in this field such as modelling, fault ana- 
lysis, etc. of six phase systems has been taken up by S.N. 
Tewari et al [8— ll]. Uma Pal et al [24] have reported on 
feasibility of twelve phase systems and their fault analysis. 
The symmetrical component transformation for six phase systems 
has been derived by Bhatt et al [18] as a direct extension of 
the three phase technique. L.P* Singh et al [22,23] have 
derived power invariant symmetrical component transformations 
for six and higher phase order systems using group theoretic 
techniques. L. Roy [20] has given an algorithm for developing 
the admittance matrix of polyphase systems. This is utilized 
to carry out the load flow and fault analysis of multiphase 
svs terns in nhasor coordinates. 
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All these studies on multiphase systems have encouraged 
electrical engineers and scientists to pursue the topic in 
more detail in order to establish the superiority of multi- 
phase systems, especially the six phase systems, over the 
conventional three phase systems. 

In the present work, the detailed representation of three 
and multiphase (six and twelve-phase) systems is done in phasor 
coordinates. This helps in taking into account various types 
of unbalances existing in the system which are, otherwise, 
left untouched. The advantage of such a representation is, 
that, it retains the physical identities of the different 
elements present in the systems. Models to analyse composite 
three-phase and six-phase systems and composite three-phase and 
twelve-phase systems have been derived. The representation is 
used to carry out the balanced as well as the unbalanced load 
flow analysis of three and six T phase systems. We now give the 
chapterwise description of the work carried out in this thesis. 

1.2 CHAPTERWISE DESCRIPTION 

Chapter 2 starts with the emphasis on representing the 
three-phase system in phasor coordinates. The detailed modell- 
ing of three-rphase elements , such as machines, different types 
of transformers, transmission lines has been given using the- 
symmetrical lattice equivalent circuits. This modelling is 
employed to carry out the balanced as well a s unbalanced load 
flow analysis of a sample power system. 
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Chapter 3 initially presents an overview of the 
possibility of multiphase transmission systems from steady 
state and overvoltages and insulation requirements point of 
view. This is followed by mathematical modelling of the 
different components such as synchronous machines, trans- 
mission lines and transformers of six. phase systems. The 
transmission line has been modelled based upon its phase 
impedance matrix, ABCD parameters and it-equivalent circuit. 
Three-phase equivalent representation of a six-rphase system 
and vice versa has been derived. This equivalent representa- 
tion is utilised to carry out the load flow studies of a sample 
power system on equivalent three phase basis. The chapter 
concludes with the load flow analysis of a composite three and 
six-phase power system. 

Chapter 4 starts with the development of detailed mathe- 
matical model of machines, transformers and transmission lines 
of twelve-phase systems in phasor coordinates. Phase impedance 
matrix representation, ABCD parameter representation and the 
it-equivalent representation of the twelve phase transmission 
line has been given. The three-phase equivalent representation 
of a twelve-phase system and vice versa have also been developed. 

Chapter 5 summarises the work reported in this thesis and 
the conclusions drawn from the results and findings of the work. 



CHAPTER 2 


THREE PHASE POWER SYSTEM: MODELLING AND ANALYSIS 

2.1 INTRODUCTION 

Traditionally the analysis of power systems is done based 
on the assumption of a balanced system. This is done by 
representing the system by a one line diagram. Although such 
an analysis gives sufficiently accurate results, it fails to 
take into account the inherent unbalances present in the system. 
There are a number of situations in practical power systems 
like very long untransposed transmission lines, bundle condu- 
ctors, large single phase loads, one, two or four phase supply, 
single pole switching, etc. which make the system unbalanced 
and warrant the need for a more detailed representation of the 
sys tem. 

From the analysis of a balanced system, various parameters 
namely, voltage magnitude, its angle, current, real and reactive 
powers at a bus are found for one phase only and those for the 
other two phases are obtained from the concepts of a balanced 
system. However, in an unbalanced system, such a procedure 
cannot be employed because of the obvious reasons . 

With a view tb analysing the power systems taking into y 
account the unbalances present in it, it becomes imperative tb 
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represent the system in detail. This requires a critical 
modelling of the various elements of a power system such as 
machines, transformers, transmission lines and loads, etc. 

This chapter deals exclusively with the modelling of three 
phase power systems. The chapter starts with the derivation 
of the final equation of a machine representing all the un- 
balances in it. After this the representation of different 
types of transformers are given based on their symmetrical 
lattice equivalent circuits. Finally the modelling of trans- 
mission lines and loads are done. This modelling of the system 
is employed to carry out the three phase load flow studies of a 
saipple system in phasor coordinates. 

2.2 MODELLING OF THREE PHASE MACHINES 

The general three phase element can be represented as shown 
in Fig. 2.1. Fig. 2.1 represents a wye-connected machine if 
E a ,Ej^,E c are a balanced set of internal voltages and nodes 
4,5 and 6 are joined to form a star point [15]« Hence, a three 
phase synchronous generator is represented by a balanced set of 
internal voltages behind proper reactances. Fig. 2,2 shows the 
self and mutual admittances for any node a in detail. 

The KCL at any node a is written as 
[(V 4 + V- V lj Y aa + <- Y ab> + <~ Y ac> 

+ ^V E a>"< V E b>5 Y ab + {<V E a>-'V E c>5 Y ac * X 4 ^ 
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or. 


■*“ Y aa Y ab Y ac Y aa ” Y ab 





V 456 




( 2 . 2 ) 


where 


V 123 = [V, V 2 Vg] 

V 456 = C V 4 V 5 V 6^ 
E 


*abc 


= [E a E b *y 


T 

T 

T 


Similarly, from Fig. 2.2, we can write equations at nodes b and 
c also. Combining all these equations, we obtain the expression 
of the following form : 


cy y 

aa ab 

Y ac Y aa 

— Y — Y V 

ab ac aa 

~ Y ab “ Y a<T 


V 123 

Y ba- Y bb 

Y bc " Y ba 

Y bb “^bc “ Y ba 

Y bb “ Y bc 


V 456 

Y ca Y cb 

- Y cc -Y ca 

“ Y cb Y cc “ Y ca 

~ Y cb Y cc^ 


^abc 



(2.3) 


Likewise, Fig. 2.3 shows the admittances at node 1 • Thus 

KCL at node 1 is : 




FIG. 2.1 SCHEMATIC REPRESENTATION OF A THREE PHASE ELEMENT 






^ 2 

v 3_„ 


2.2 SCHEMATIC REPRESENTATION OF SELF AND MUTUAL 
ADMITTANCES AT ANY NODE A OF A THREE PHASE 
MACHINE. 
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FIG. 2.3 SCHEMATIC REPRESENTATION OF SELF AND MUTUAL 
ADMITTANCES AT ANY NODE 1 OF A 3-0 MACHINE 


la 



FIG. 2.4 THE INJECTED CURRENT AT THE NEUTRAL OF A 
THREE PHASE GENERATOR. 
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•^ 1 -<V E a ) } Y aa + [ v l-<V E b>j (- Y ab> + ^ Y ac> 


+ < V l- V 2) Y ab + ' V l- V 3> Y 


ac 


= I. 


(2.4) 


or 


^aa ~^ab ~^ac “^aa Y ab Y ac *“^aa ^ab ^acl 


123 

456 


'abc 


- Il 
(2.5) 


Proceeding in the same fashion, we can form equations similar 
to eqn. (2.3) by writing KCL at nodes 2 and 3. Combining these 
three equations results in the following equation : 


aa 

“ Y ab 

- Y ac 

- Y aa 

Y ab 


ba 

Y bb 

Jk 

cr 

o 

Y ba 

- Y bb 

Y bc 

ca 

* Y cb 

Y oc 

Y ca 

Y cb 

- Y cc 


’ Y aa Y ab Y ac 


V 123~ 


IT 

I 1 

Y ba - Y bb Y bc 


V 456 

ss 

I 2 

Y ca Y cb- Y cc 

1 

J^abc 


V 


(2.6) 


Combining eqns. (2.3) 

and (2.6) 

, we 

get 









J 

Y abc ~^abc 

”^abc 


V 123 


I 123 


~^abc ^abc 

^abc 


V 456 

C 


I 456 



— i 


^abc^ 

, 

— - 


where, 


Y aa 

- Y ab 

~ Y ac 


Y . - 

abc 

“^ba 

Y bb 

" Y bc 


i 

”^ca 

- Y cb 

Y cc 


(2.7) 


( 2 . 8 ) 
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and 


*123 ~ ^1 X 2 I 3^ T » 
I 456 = ^■ I 4 I 5 I 6^ T 


Eqn. (2.7) is the concise form of nodal relationships for the 
current and voltage reference directions shown in Fig. 2.1. 

If the parameters of machine are given in the component 
form, the phase admittance matrix Y a};)C is obtained by 

■ T 


Y abc “ T s Y °’ 1 ’ 2 T s 


(2.9) 


when 


T s = 


f3 


1 1 
1 


o c 2 a 


a a 


a = 3fl = 1/120 ° * e j2it/3 = _o,5+j 0 .866 
and 

Y 0flf2 = diag. [y o ,y 1 ,y 2 ] 

After using eqn. (2.9), we obtain 


abc 


_ 1 


( Y 0 +y l +y 2 ) ( y o +ay l +a2y 2 > ^ y o +a2y l +ay 2 > 

(y Q +a Yi+ay 2 ) (y 0 +y i +y 2^ ( y o +0Cy i +a y 2^ 

( y 0 +ayi+a 2 y 2 ) (y 0 +a 2 Y;j+aY 2 ) (y 0 +y i +y 2^ 


( 2 . 10 ) 
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where Y 0 »y^ and are the zero, positive and negative 
sequence admittances of the synchronous generator respectively. 
For balanced machine interval voltages, 

E b - cc 2 E a and E c = a E a 

If in Fig. 2.1, nodes 4,5 and 6 are joined to form a neutral 
point N, then we have 

1 n =1 4 + h + h (2- 11 ) 

and 

V N = V 4 = V 5 - v 6 < 2 - 12 > 

where 1^ is the injected current at the neutral N and is 
the neutral voltage. 

If S a ,Sk and S c denote the complex powers in phases 
a,b and c respectively, then 


V E a> V N +E b and V E c 


are the interval terminal 


voltages and 


I- = 


<V E a> 


> Ih = 


(v N + E b )' 


= 


(¥ E C > 


(2.13) 


Substituting the value of ^abc from eqn. (2.10) into eqn. 
(2.7) and also the values of E b and E c and taking 

I 4 = 1 a' *5 = I b > I 6 = V 
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eqn. (2.7), after collection of common terms for and E , 
reduces to the following : 

(2.14) 


From the eqn. (2.11) , we note that an additional eqn. can be 
written for node N. Hence, after proper substitution from 
eqn. (2.10) in eqn. (2 Jfa) and adding last three equations of 
(2.14) gives : 



~ y o 

~ y l 

Y abc 

- y o 

-a 2 yi 


~ y o 

“ ay l 


y o 

y l- 

~ Y abc 

y o 

2 

a- y x 


y o 

ay l 


v i 


1 

• 1 

V 2 


I 2 

v 3 


I 3 

V N 


.' 4 . 

S a/‘V E a) 

E a 


S b^ V N +E b^ 

_ J 

1 

_fc/(V E c)_ 


[-Y, 


•Y« “ 


y o 3y o 03 


V, 


V, 


V, 


N 




(2.15) 


* 


* 




Noting that E K = <xE and = a E - the last three right hand 
jo a c a 

Side terms of eqn. (2.14) become 




i I j.-, 


& +E a 


* K 

V N +aE a 


and 


* 2 * 
V“ E a 


respective ly. 


.ifjross multiplying the last three rows of eqn. (2.14), we obtain 


efi 

41 
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T [-(Y 0 + Yl +y 2 > V 1 (V*+E*)-(y 0 +ay 1 + a 2 y 2 )V 2 (V*+E*) 




*r 




-ly 0+ a yi +ay 2 )V 3 (V N +E a )]+ Y 0 V N (V N+ E a ) +Yl E a (V N+ E a ) = S a 

(2.17) 

| [-(yo+^Yi+ayz) V 1 (V*+aE*)-(y 0 +y 1 +y 2 )V 2 (V*+aE a ) 


-(y 0+a y 1+ a-y 2 )V3(V* + c t Ep] + y 0 V N (V^o ! E;) + o ! ' : y 1 E a (v;-H I E;) = S’ b 

(2.18) 


■fr 




% [ -( y o +ay 1 +a 2 y 2 ) V x ( v*+a 2 ^) -( y 0 +a z y 1 +ay 2 ) V 2 ( V~+a 2 E~ ) 


r* 2„*- 


“ ( Y o + Y 1 + Y2 ) ^3 ( Y o V N ^ ^N +<x ^ a ) +ot Y 1^ a ^ ^N +a ^ a ^ = S c 

(2.19) 


Adding eqns. (2.17), (2,18) and (2.19) and simplifying we 
obtain 


-yo v i v N-Vo v 2 v N-y 0 v 3 v N +3 yo v N v N-yi v i< 

-ayi v 2 E a -o 2 yiV 3 E a+ 3yiE a E a = S* + S* + S* (2.20) 

The first four terms of eqn . (2.20) are times 1^. Hence, 
under usual conditions, 1^=0 and eqn. (2.20) reduces to 
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ry 

[-y x -ay 1 -a y x 3 y-^ 


V J 

V, 


V, 


* J* * 

S +S.+S 

3 D C 


( 2 . 21 ) 


Suppose the generator is earthed through an admittance y^Q as 
shown in Fig. 2.4, then 


I H ~ I 4 +I 5 +I 6 +y N0 V N 


= [-y 0 -y c ~y 0 3y n +y, 


o^KO 


0] 


V ] 

V. 

z 

V. 


V, 


N 


a Ji 


( 2 . 22 ) 


(2.23) 


Substituting the value of in eqn. (2.23) and putting 

1^ = 0, we get 

lv i, 

|2 


i-Y 1 -ayx -a 2 y x 3yJ 


V, 


w s . v NO fv N r 


(2.24) 


The second term on the R.H.S. of the eqn. (2.24) is negligible 
as compared to the first, term. Hence, the final machine 
equation, with 1^ = 0, is the combination of the first three 
eqns. of (2.14), eqn. (2.23) and eqn, (2.24) : 
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-Vl ~ Y o 

. 

. V 1 



Y 2 

Y abc Vi -Y 0 

’ 

V 2 


*2 

‘“Vi -Y 0 


V 3 

= 

*3 

~Y 1 “ a Y 1 -a 2 Y 1 3y x 0 


E a j 


_ * j * 

IS /E a 

Jo ~ Y o -V 0 3 Yo + Yno 


V N 


0 

1 



— ^ 


rmm 


; 


(2.25) 


It may be noted that the machine model eqn. (2.25) makes use of 
the total power and not the individual phase powers. The 
injected nodal currents 1^,12,13 are zero for the synchronous 
generator if there is no local load. However, if there is a 
local load, they are not zero and one third of the' total local 
load is assigned to each of the phases. 


2.3 MODELLING OF THREE PHASE TRANSFORMERS 


In this section, we shall first derive the symmetrical 
lattice equivalent circuit of a single phase transformer with 
off-nominal tappings on primary as well as on secondary sides. 
This lattice equivalent circuit will then be employed to 
represent different types of three phase transformers [7]. 

2.3.1 Symmetrical Lattice Equivalent of one-phase Transformers 

Fig. 2.5 illustrates a single phase representation of a 
transformer in per unit form as an ideal transformer of turns 



FIG. 2.5 SCHEMATIC REPRESENTATION 
PHASE TRANSFORMER. 


OF AN IDEAL SINGLE 


y/d*t) 



FIG. 2.6 SYMMETRICAL LATTICE EQUIVALENT 
FIG. 2.5. 


CIRCUIT OF 
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ratio (l+t):l in series with an equivalent leakage admittance 
y. From Fig. (2.5) and the relationships in an ideal transfor- 
mer, we can write 


V o = (l+t)V* 

o 

i r =-d+t)i s 
The terminal relationships are 


= y(V r -V' ) = yV r - 


1+t 


i 


s 



d+t ) 2 


yv r 

TT+t) 


(2.26) 


(2.27) 


If we do not assume any earth connections and represent the 
nodes of the transformer by j,k,p,q, then 


V s - Vj-Vk and v r = V p- V q ( - 2 ' 28 '> 

The injected currents into each node for an n-node network are 
expressed in terms of these nodal voltages by equations of the 




E ( V.-V )y . + 

m J m jm (j_ + t) 


5 (VV-nitT ( VV 

(2.29) 


form 
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where the summation is over the set of all nodes m connected 
to node j excluding the set k,p,q. Expanding the eqn. (2.29) 
for each node of transformer, we obtain 




(1+t) 2 k 


y ,, 

Tl+tJ V P 


+ TT+tJ V q + l ( - y jm )V m 


*k " 


(1+t) 


Y 5 v i + ( 5 y km + 7rrT5 )v k + n+tT V P 


m 


(1+t)' 


~ TT+TJ v q + s ( -Ykm )v m 


(2.30) 


(2.31) 


i p=- tmr v i + iw \ + v P m + y>V y V z( - y p» )v »n . 

(2.32) 

1 q = i4 v j - imy v k - yV P + y qm +y) V * ( -v )v m 

, (2.33) 

The nodal admittance relationships as seen from eqn. (2.30) to 
eqn. (2.33) may be written in the form 

* ' 

I = YV + Y* V (2.34) 


The equivalent circuit shown in Fig. 2.6 is deduced from the 
algorithm of forming matrix Y (by inspection). Thus it 
represents a single phase tapped transformer. The voltage and 


21 


y/of n 



FIG. 2.7 GENERAL SYMMETRICAL LATTICE EQUIVALENT CIRCUIT 
OF A SINGLE PHASE TRANSFORMER. 



FIG. 2.8 SYMMETRICAL LATTICE EQUIVALENT CIRCUIT OF A 

3-0 STAR-STAR TRANSFORMER AS A LINEAR GRAPH. 
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current relationships between the nodes of this symmetrical 
lattice equivalent circuit are same as those in Fig. 2.5* 

From Fig. 2.6 we can easily obtain the circuit shown in 
Fig. 2.7 [7] which represents the general symmetrical lattice 
equivalent circuit of a single phase transformer, where both 
primary and secondary windings may have either actual or 
equivalent variable turns ratio a and (3 or both. 

2.3.2 Representation of Three-Phase Transformers 

The general symmetrical lattice equivalent circuit derived 
in the Section 2.3.1 is employed to represent the equivalent 
circuit of different types of three phase transformers, namely, 
star-star, delta-delta and star-delta in this section. 

a) Star-Star Transformer : 


The lattice equivalent circuit of a three phase star-star 
transformer as a linear graph is shown in Fig. 2.8 in which 
parallel transformer windings are taken to represent the 
equivalent single phase transformers. The circuit is drawn in 
analogy with that of Fig. 2.7 with taps on only one (primary) 
side since in practice either a or p would be 1 p.u. The 
’ equivalent circuit may be described by the connection Table 2.1. 



b) Delta-Delta Transformer s 

The lattice equivalent circuit of the three phase delta- 
delta transformer as a linear graph with the same convention 
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Table 2.1 

Connection table for Y-Y transformer with 
a = l+t a , p = 1 p.u. 


Admittances 

Between Nodes 

y/a 2 

N-A, N-B, N-C 

Y 

n-a, n-b, n-c 

y/a 

A-a, B-b, C-c 

i 

. i 

-y/a 

j 

n-A, n-B, n-C; N-a, N-b, N-c 

3y/a 

N-n 

< 


that parallel windings may be considered to represent single 
phase transformer is constructed as in Fig. 2.9. The conne- 
ction table 2.2 describes the admittances between different 
nodes for such a transformer with the assumption of off- 
nominal tappings only on the primary side. 

c) Star-delta Transformer : 

Using the same conventions and techniques as above, the 
lattice equivalent circuit of a three-phase star-delta 
transformer as a linear graph is shown in Fig. 2.10. The 
convention used for numbering the nodes and thus identifying 
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Table 2.2 

Connection table for a delta-delta transformer 
with a = f3(l+t a ),£ = Y3 p.u. 


Admittances 

Between Nodes 

y/a 2 

A-B, B-C, C-A 

y/3 

a— b , b— c , c— a 

2y/V*3a 

A-a, B-b, C-c 

-y/V3a 

A-b, B-c, C-a; a-B, b-C, c-A 


opposite sides of the symmetrical lattice networks is as 
follows : 

% 

A-N/c-b; B-N/a-c; C-N/b-a 

The connection Table 2.3 shows the admittances between various 
nodes. Also, for taps only on primary side, t^ = 0. 

2.4 MODELLING OF THREE-PHASE TRANSMISSION LINES 

A* «*transposed three phase transmission line is described 
by its phase impedance matrix for a three phase order as 
follows : 





FIG. 2.9 SYMMETRICAL LATTICE EQUIVALENT CIRCUIT OF 
A 3-0 DELTA- DELTA TRANSFORMER AS A 
LINEAR GRAPH. 


B a 



FIG. 2.10 SYMMETRICAL LATTICE EQUIVLENT CIRCUIT OF 
A 3-0 STAR-DELTA TRANSFORMER AS A 
LINEAR GRAPH. 
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Table 2.3 

Connection table for a star-delta transformer with 
a = (l+t a ) , § = lf3 p.u. 


Admittances 

Between Nodes 

y/a 2 

A-N, B-N, C-N 

y/3 

a-b, b-c, c-a 

y/f3a 

: 

A-c , B-a, C-b 

-y/f3a 

A-b, B-c, C-a 



‘s 

Z ml 

Z m2 

‘m2 

Z s 

Z ml 

ml 

\q. 

Z s^ 


(2.35) 


3 

This coefficient matrix Zp can be diagonalized by the symmetri- 
cal component transformations. 


Although, short lines can be represented by eqn. (2.35), 
the adequate representation of the medium and the long lines 
is obtained by nominal and equivalent n circuit representations. 
Fig. 2.11 shows schematically the series admittances Y a k C of a 
3-phase transmission line between buses 1,2,3 and 4,5,6 and 
self and mutual admittances at any node 1. 
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Applying the KCL at node 1 , we get 


' V l- V 4> Y aa + ( V l- V 5H“Y ab )+(V 1 -V 6 )(-Y ac ) + (V 1 -V 2 )Y ab+ (V 1 -V 3 )Y ac =I 1 


fv -Y ~Y -Y Y Y ""I 
aa ab ac aa ab acj 


= [In 3 


(2.36) 


(2.37) 


Similarly, writing the KCL at all the nodes and combining them 
we get the following result 



(2.38) 


(2.39) 


where Y abc is given by eqn. (2.8). 
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For equivalent % representation of the transmission 
line, we have also got to consider the shunt parameters at 
each end of the line as shown in Fig. 2.12. Hence, the final 
nodal voltages and currents relationships for the considered 
transmission lines is given by [7] 



(2.41) 


If we have got two parallel lines, the couplings between the 
two can be considered as shown in Fig. 2.13. Each element in 
Fig. 2.13 represents a (3x3) admittance matrix. The bus 
admittance matrix for these two lines is given by 




FIG. 2.11 SCHEMATIC REPRESENTATION OF SELF AND MUTUAL 
ADMITTANCES OF A 3-0 TRANSMISSION LINE AT A 
NODE. 


Yab/2 



FIG. 2.12 THE CAPACITIVE SUSCEPTANCES AT ONE END 
OF A 3-0 TRANSMISSION LINE. 
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k 

1 

m 

n 

k 

n r~r * -M 

Ykl + Ykl,g 

“ y kl 

y kl ,mn +y kl,mn,g 

~ y kl,mr 

1 

-y k i 

Ykl + Ykl,g 

~ y kl ,mn 

y k 1 , mn + ; 





y kl ,mn, g 

m 

y kl,mn +Y kl^mn,g 

“ y kl,mn 

y mn +y mn , g 

-y 

7 mn 

n 

~ y kl ,mn 

y k 1 , mn +y kl , mn , g 

- y mn 

y mn +y mn , c 


(2.42) 


2.5 THREE PHASE LOADS 


If three phase static loads are present in a power system, 
they are assumed to be decoupled and the admittance equivalent to 
one third of the load is added to the corresponding diagonal 
element of the system admittance matrix. It is calculated from 
the complex power and prefault nodal voltage. The load admi- 
ttance is obtained from the following eqn. [20] 


Yk = 


Pk-j^k 




(2.43) 


where k is any node. 

Eqn. (2.43) is solved for each phase which becomes a bus in 
this case using the phase power and the phase nodal voltage. In 
the absence of any voltages, the three phase balanced voltage 
with unity magnitude may be used. 
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In case of a single phase load, its equivalent admittance 
is added to the corresponding diagonal element just like a 
shunt appearing between one of the nodes and the ground. The 
compensating reactors are treated like static load. 

2.6 SAMPLE NETWORK 

The single phase and three phase (balanced as well as 
unbalanced) load flow analysis of the sample power system L ,c Vl 5 _2 
network as shown in Fig. 2.14 was carried out. The bus 
numbering sequences for single phase and three phase load 
flow have been shown separately in the figure. The necessary 
data for carrying out the analysis of the sample network are 
given in Table 2,4. The familiar Gauss-Seidel iterative 
technique was used in the study. The results obtained from 
the single phase load flow analysis with a tolerance of 
.00001 are shown in Tables 2.5 - 2.6. The computations were 
carried out on the DEC-1090 System at I.I.T., Kanpur, India. 

Balanced Three Phase Load Flow Analysis : 

The three phase balanced load flow analysis of the sample 
network shown in Fig. 2.14 was done in the phase frame of 
reference. The procedure for assembling the polyphase network 
nodal admittance matrix given by Roy [20] was employed. It 
may be noted, that, in three phase load flow studies, every 
phase at each busbar including the neutral of generator 
becomes a node. The main steps for carrying out the load flow 
analysis are : 
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i) Assign each phase of each busbar including neutral of 
the generator a node number. 

ii) Form the system multiphase nodal admittance matrix using 
the representation of elements in the phase frame of 
reference . 

iii) Form and then solve the nodal performance equation, 

I = YV, using some iterative technique (like Gauss- 
Seidel iterative technique in the present case) 

iv) Compute the voltage magnitude, its angle at each bus 
and line flows, etc. 

The results obtained from the balanced three phase load 

flow analysis, whose flow chart is given in Appendix , are 

given in Tables 2.7 - 2.8. 

Unbalanced Load Flow Analysis : 

The following unbalances were simulated in the sample 

network to evaluate the performance of the system. 

i) The loads at node number 13,14,15 were unbalanced to 
110 percent, 100 percent and 90 percent of their base 
case values. 

ii) The delta winding of the transformer connected between 
the nodes 4,5,6 and 16,17,18 was converted to an open 
delta. 

iii) The off-nominal tappings of the secondary side of the 
transformer connected between nodes 4,5,6 and 16,17,18 was 
set at -0.020 p.u. 
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The results obtained from the load flow analysis after 
making the changes as above are shown in Tables 2,7 - 2,8* 

2.7 CONCLUSIONS 

In this chapter, the necessity of representing the three- 
phase system in phase frame of reference has been emphasised. 
Consequently, the detailed mathematical models of the three 
phase generator, transformer, transmission line and load has 
been done. It has been shown that, more realistic performance 
of the system can be studied by doing phasor coordinate analy- 
sis of the system. 

The following conclusions can be drawn, from the load flow 
analysis of the sample network done in the previous section : 

i) The results obtained from the single phase load flow 
analysis and the balanced three phase load flow analysis 

of the sample system are alike. This proves the validity 
of the single phase load flow analysis based on one line 
diagram when the system is balanced. 

ii) When the system becomes unbalanced, single phase load flow 
fails to give accurate results as seen from the unbalanced 
load flow analysis. 

iii) The unbalances in the magnitudes of voltages of different 
phases at a bus are considerable. 

iv) It is seen, that, the line flows in the three phases 
between two buses are quite different from each other 
during unbalanced conditions. 
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The above results and the conclusions establish the 
need for conducting the load flow analysis of the system in 
phasor coordinates, especially, when there are unbalances in 
the system. 
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(4-6) (31-33) 



16.. 2.14 THE SAMPLE NETWORK WITH BUS NUMBERING SEQUENCE FOR 
SINGLE PHASE AND THREE PHASE LOAD FLOW STUDY. 



37 


Table 2.4 


Circuit characteristics for sample system of Fig. 2.14 
on 100 MVA base/phase and rated KV 


From To 
Bus Bus 


+ Sequence 
R 


Zero sequence 


X. 


1) Generation 


U-3) 

(28-30) 

0.0 

0.0967 

0.0 

0.0467 


(4-6) 

(31-33) 

0,0 

0.17 

0.0 

0,0850 

— 

(7-9) 

(34-36) 

0.0 

0.17 

0.0 

0.0850 

- 


2) Single circuit 3-0 Lines, : 


(25-27) (22-24) 
(25-27) (13-15) 
(16-18) (13-15) 
(25-27) (16-18) 
(19-21) (16-18) 


0.05 0.20 
0.10 0.40 
0.15 0.60 
0.10 0.40 
0.05 0.20 


0.033 

0.066 

0.100 

0.066 

0.033 


0.155 

0.310 

0.465 

0.310 

0.155 


'<) Double circuit 3-0 Lines 


10-12 

10-12 

22-24 


'13-15) 0.0375 0.150 0.00625 
25-27) 0.01875 0.075 0.00313 
,19-21) 0.025 0.100 0.00415 


0.11625 

0.05813 

0.0775 


!■) Transformers : 

(1-3) (10-12) 

(4-6) (16*18) 

(7-9) (22-24) 


Type of connection 

& - Yx 
& - 'h. 

& - 'H. 


0.580 

1.160 

1.740 

1.160 

0.580 


0.435 

0.2175 

0.2900 


0.0193 

0.0387 

0.0580 

0.0387 

0.0193 


0.0036: 

0.0018C 

0.0024C 


Reactance X 

0.0533 

0.12 

0.16 


Table 2.5 

Voltage solution for single phase load 
flow of the network of Figure 2.14 


Bus Number 

Voltage Magnitude 
in p.u, 

Angle in 
degree 

1 

1.043 

-2.745 

2 

1.027 

-2.642 

3 

1.005 

-4.132 

4 

1.030 

-1.465 

5 

1.023 

-2.967 

6 

1.022 

-3.746 

7 

0.989 

-1.265 

8 

0.996 

-4.943 

9 

1.018 

-2.089 

10 

1.070 

0.4TO» 

.11 

1.040 

-1.378 

12 

1.020 

-0.986 


Table 2.6 


Power Flows 


From Bus 

To Bus 

Real Power 
(MW) 

Reactive Power 
(MVAR) 

10 

1 

276.554 

98.746 

11 

2 

179.295 

34.193 

12 

3 

139.836 

30.574 

8 

9 

-28.761 

-30.792 

5 

9 

-0.498 

-7.864 

5 

6 

+14.231 

49.243 

6 

9 

-22.902 

-12.962 

6 

7 

-93.807 

-24.024 

5 

4 

-90.287 

-9.356 

9 

4 

-185.176 

-52.262 

7 

8 

+88 .044 

+22.575 






Table 2.7 

Voltage solution for three phase load flows of the 
network of Fig. 2.14 


Bus No. 


Balanced Case 


Unbalanced Case 


( lvl 1 
(p.U.) 

6 

(deg. 

1.043 

—4.616 

1.043 

-124.686 

1.043 

115.446 

1.026 

-2.831 

1.026 

-122.888 

1.026 

117.189 

1.005 

-5.617 

1.006 

-125.675 

1.006 

114.401 

1.031 

82.696 

1.031 

-37.231 

1.031 

-157.281 

1.023 

81.483 

1.023 

-38.444 

1.024 

-158.489 

1.021 

83.190 

1.021 

-36.740 

1.022 

-155.768 

0.988 

79.896 

0.988 

-40.035 

0.999 

-160.060 

0.996 

80.687 

0.996 

-39.242 

0.996 

-159.270 

1.018 

81.513 

1.019 

-38.415 

1.019 

-158.458 

1.070 

0.000 

1.070 

-120.072 

1.070 

120.072 

1.040 

2.621 

1.041 

, -117.436 

1.041 

122.629 


i\ 

8 

8!SfS 


121.907 
118.158 
0.0 


0.000 

0.000 


T5* 

• < 

c — 

• 

Q 

(deg.) 

1.040 

-4.951 

1.034 

-124.585 

1.044 

115.604 

1.029 

89.329 

1.014 

-29.109 

1.040 

-149.156 

1.002 

-6.033 

0.996 

-125.508 

f 1.007 

114.587 

1.023 

83.075 

. 1.034 

-37.435 

1.018 

-157.618 

1.013 

81.764 

1.029 

-38.707 

1.007 

-158.672 

1.011 

84.550 

1.031 

-37.209 

0.981 

-157.437 

0.981 

80.615 

0.994 

-40.263 

0.963 

-160.759 

0.988 

81.307 

1.000 

-39.487 

0.976 

-159.857 

1.011 

81.959 

1.023 

-38.620 

1.002 

-158.874 

1.070 

0.000 

1.070 

-120.072 

1.070 

120.072 

1.040 

95.850 

1.041 

-24.207 

1.041 

-144.215 

1.020 

-1.862 

1.021 

-121.919 

1.021 

118 . 146 

1.021 

0.000 

0.000 

0.000 

0.000 

0.000 

, 






Table 2.8 


Power Flows 


> From Balanced Case 

is Bus 


Real Power Reactive Power 
(MW) (MVAR) 


28 

92.919 

' 33.974 

29 

92.901 

33.691 

30 

93.129 

33.860 

31 

59.656 

11.698 

32 

59.761 

11.889 

33 

59.565 

11.802 

34 

39.667 

10.098 

35 

39.748 

10.225 

36 

39.586 

10.192 

22 

-9.507 

-10.402 

23 

-9.557 

-10.548 

24 

-9.412 

-10.523 

13 

0.159 

-2.285 

14 

0.162 

-2.278 

15 

0.154 

-2.281 

13 

-4 . 79 3 

-3.654 

14 

-4.771 

- 3.597 

15 

-4.833 

-3.610 

16 

7.367 

-4.422 

17 

7.337 

-4.508 

18 

7.423 

-4.483 

16 

31.465 

8.041 

17 

31.452 

8.031 

18 

31.470 

8.031 

13 

-30.368 

-3.061 

14 

-30.391 

-3.122 

15 

-30.320 

-3.109 

25 

-62.132 

-17.570 

26 

-62.207 

-17.786 

27 

-61.981 

- 17.742 

19 

-29.050 

- 7.707 

20 

-29.055 

-7.711 

21 

-29.036 

-7.712 


I 

V 

L. 

\ 

% 

? 

If 
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Unbalanced Case 


Real Power 
(MW) 

Reactive Power 
(MVAR) 

99.341 

38.021 

90.083 

43.307 

90.026 

32.727 

70.727 

10.923 

53.109 

18.730 

55.140 

-0.499 

43.753 

12.556 

37.685 

16.187 

37.567 

9.072 

-8.175 

-11.052 

-10.042 

-9.972 

-10.526 

-11.478 

-0.331 

-2.374 

0.180 

-2.311 

0.641 

-2.087 

-7.544 

-2.209 

-4.036 

-4.757 

-2.786 

-0.128 

11.106 

-6.198 

6.234 

-2.962 

4.827 

-9.080 

36.144 

5.818 

30.231 

10.061 

27.867 

1.061 

-30.622 

-5.317 

-31.143 

-1.935 

-29.355 

-5.885 

- 57*591 

-20.641 

-63.763 

-15.477 

-65.209 

-23.952 

-24.521 

-10.560 

-30.296 

-5.50 

-32.542 

-14.454 






CHAPTER 3 


SIX PHASE POWER SYSTEM: MODELLING AND ANALYSIS 
3.1 INTRODUCTION 

t 

Three phase power systems are the most widely used in the 
present day electric power systems* Although high phase order 
machine construction and power transmission have been under 
active consideration, specially during the last 10 years, still 
three phase systems continue to be in use. However, as the 
demand for power is increasing at a phenominal rate, multiphase 
systems are bound to become a reality in near future. 

All the advantages of the multiphase systems stem from 

the fact that the phase angle between two phases is less 

vo U&vws 

resulting in lower adjacent phase awptees and also, lower line 
to line voltages between two adjacent lines. 

Out of the higher phase orders proposed so far in papers 
and conferences (6,9 and 12), the six phase systems appear to 
be most promising. To carry out the planning and operation 
studies of six phase systems, their mathematical modelling 
should be known beforehand. Mathematical modelling is an easy 
and flexible means to understand the behaviour of the system. 

This chapter deals mainly with the modelling aspect of 
various elements present in a six phase system. The chapter 
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starts with an overview of the feasibility of multiphase 
systems with special reference to six phase systems. The 
modelling starts with the development of an equation for a six 
phase machine including all the unbalances present in it. 

After that, representations of three/six-phase and six phase 
transformers have been proposed based upon their symmetrical 
lattice equivalent circuits. In the next section, the repre- 
sentation of a six-phase transmission line based upon its 
phase impedance matrix, rc equivalent circuit and ADCD para- 
meters has been given. To carry out the analysis of a compo- 
site three/six phase system, three phase equivalent represen- 
tation of a six phase line and vice versa has been derived. 

This is followed by the modelling of six phase loads. Finally, 
the chapter concludes with the load flow analysis of a sample 
system along with the results and discussions. 

3.2 FEASIBILITY OF MULTIPHASE SYSTEMS 

The feasibility studies of multiphase systems, in parti- 
cular, six phase systems, are available in the literature 
[2-6,8-11] with regard to their steady-state operation, over- 
voltages, insulation requirements, power flow, reliability, 
etc. It has been found [2-3] that in several features, six 
phase systems possess better characteristics than that of the 
three phase systems. The summary of the relative merits and 
demerits [2-3] of the multiphase systems is given below. 
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(a) Power Transfer Capability ; 

More power can be transferred with multiphase systems 
and smaller line and tower dimensions are required. It has 
been found out that, if an existing double circuit three phase 
line is converted to a six phase line with the same adjacent 
phase-phase voltages in both the cases, the power transfer 
capability gets increased by 73.2 percent. 

(b) Conductor Clearance : 

If the system voltage is defined with phase— ground voltage 
as references, we find that adjacent phase-phase voltage de- 
creases with phase order. Hence, line conductor spacing 
reduces with increased phase order.. But there is a limit to 
reduction of spacings because of motion of individual condu- 
ctors owing to ice, wind, faults, etc. 

(c) Transposition : 

It is difficult to transpose multiphase lines unlike the 
three phase lines which can be freely transposed. The only 
practical transposition for multiphase lines is obtained by 
rotating the entire conductor array in steps over the length 
of the line in a cyclic manner. 

(d) Thermal Loading : 

If the thermal loading is the criterion for circuit 
rating, the capacity increase is proportional to the number 


\ 


45 


of phases as the thermal loading is directly proportional to 
the number of phases. 

(e) Surge Impedance Loading : 

The Surge Impedance Loading of a transmission line is 
almost proportional to the phase order increase and reaches 
saturation beyond six-phase. 

(f) Current Unbalance : 

Current unbalances in multiphase lines are less than 
that of three phase lines. Hence, six phase circuits are 
better balanced compared to two three phase circuits with the 
same conductor configurations. Hence, transposition may be 
unnecessary in the former case. 

(g) Electric Fields : 

Shield wires may be required for high phase orders 
to compensate the effect of decrease in maximum surface 
electric field and increase in maximum ground electric field. 

(h) Radio and Audible Noise : 

The performance of multiphase systems is better than 
their three phase counterparts as far as radio and audible 
noise is concerned, 

(i) Fault Overvoltages : 

For a six phase system, fault overvoltages are slightly 
higher than those in a comparable three phase system. 
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(j) Switching Surges ; 

Although phase to ground switching surges are approximately 
same for three and six phase systems under similar conditions, 
phase to phase surges increase with phase order. This may 
impose an upper limit to achievable phase order. 

(k) Rate of Rise of Recovery Voltage 

For the same short circuit MVA, Rate of rise of recovery 
voltage across the breaker terminals during normal opening is 
less in higher phase order systems. 

( l) Lightning Performance : 

Due to the reduced dimensions of six phase lines, the 
lightning strokes to the line are about 20 percent less than 
those in three phase systems. However, the probability of 
phase-phase flashover is much higher on a six-phase line. 

Thus, assuming effective shielding, the overall lightning 
performance of high phase order lines is comparable to 
conventional three phase line, 

(m) Terminal Insulation Level : 

The terminal insulation level is slightly higher for six 
phase systems than for three phase systems or for systems of 
12 phase and higher. 

(n) Reliability i 

The reliability study of a six phase line [5] obtained 
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from a double circuit three-phase line based upon load flow 
analysis, stability analysis and stochastic reliability 
reveals that, 

(i) Voltage regulation and efficiency are better with six 
phase lines 

(ii) A six phase line seems to be more stable than a double 
circuit three-phase line 

(iii) More energy demand can be met at a load point by 
converting a double circuit three phase line to a six 
phase line 

(iv) The loss of a six-phase lii ay have greater conse- 
quences on power system than the loss of a double 
circuit three phase line. 

Although, multi-phase lines are not being used in 
practice, developments and tests in this field have shown 
that a multiphase transmission system is realizable and may 
find applications in narrower rights of way, uprating and 
multi-circuit configurations in entire transmission range of 
voltages. 

3.3 MATHEMATICAL MODELLING OF SIX PHASE MACHINES 

Although, the topic of feasibility of a six phase machine 
is still in its infancy, it has been proposed to develop six 
phase machines for high power applications. 
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In this section, a six phase synchronous generator 
has been modelled with a grounded reactance Y^q in the same 
manner as the three phase generator . The general six phase 
element can be represented as shown in Fig. 3.1* The equation 
of this element in the impedance form can be written as 


v p- + e p ■ z p x p' + v p 5 


(3.1) 


where 


V 6 , 

P' 


= t V 7 V 8 V 9 V 10 V 11 V 12^ ’ 


E = 
P 

r 6 


[E a E b E c E d E e E f ] T , 

T 


V p = ^ V 1 V 2 V 3 V 4 V 5 V b3 > 


Zp = 6x6 impedance matrix = [Yp]**^ , 


and 


6 r 

!pt = [I 7 I 8 I 9 I xo Z 11 I 12 -‘ 


Thus, the nodal equations for the voltage and current dire- 
ctions shown in Fig. 3.1 can either be written from the eqn. 
(3.1) or derived from the symmetrical lattice equivalent 
circuit in the similar manner to three phase as follows : 


"if: 


Y^ -Y?j 

P : 


p p p 

; / ' 

I 6 i 


^6 ' 

P' ■ 


P P P 

mum 


where, I* - [l x I 2 I 3 I4 I5 I 6 1 T 

r i 


V 6 

P 

v 6< 



( 3 . 2 ) 


j* 
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Fig. 3.1 can represent a six phase machine if H a , E^,E c# E c j, 

E g and E^ are a balanced set of internal voltages and nodes 
7-12 are joined to form a star point. Hence, a six phase 
generator is represented by a balanced set of internal 
voltages behind proper reactances. For a balanced set of 
internal voltages, 

„ * * 

E, = oc E , E ~ — ccE • E. — — E , E = — <x E , Ex = &E 
b a c a* d a' e a f a 

i .e. , 

Ep ■ [l a* -a -1 -a* a] T E a (3.3) 

where, a = e 32 n :/6 = 0.50 + j0.866 = 6lTl 


If the nodes 7-12 are joined to form a star point, then 


*N ~ I 7 +I 8 +I 9 +I 10 + - I 'll +I 12 


(3.4) 


and 


'N 


V ? = V 8 = Vg = V 1Q =* V n = V 12 


(3.5) 


Hence, internal voltages become 


4 - [E a +V N . V V N • • • * W T 

£ 

and the current vector I°i is given by 


(3.6) 


* 


X a - » 5 

<W 


*b = 


<V E b> 


T 5 • •• » I. = 


<W* 


(3.7) 
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Y 6 = 
P 


s 

Y ml 

Y ra2 

Y m3 

Y m4 

Y m5 

m5 

Y 

s 

Y ml 

Y m2 

Y ra3 

Y m4 

m4 

Y m5 

Y s 

Y ml 

Y m2 

Y m3 

m3 

Y m4 

Y m5 

Y s 

Y ml 

Y m2 

m2 

Y m3 

Y m4 

Y m5 

Y 

s 

Y ml 

ml 

Y «2 

Y m3 

Y m4 

Y m5 



(3.8) 


J 


Here is obtained by 

r 


Y 6 = [T ] [Y 6 ][t*] T 

p L s J LI comp JLi s J 


(3.9) 


where T is a six phase symmetrical component transformation 
s 

matrix [21,22] given by 


T s* ~ 

5 V6 


1 

1 

1 

1 

1 

1 


* 

a ~a 


1 

-1 


* 


-a 


-1 


* 


-a 


a 


1 


-1 

1 

-1 


-a 


-a 


a 


-a 


-1 


-a 


-a 


-a 


a 


-a 


*1 


a 


,6 s 

comp 


and Y v .« diag. CY 0 tY 1 f Y 2 ^Y 3 #Y 4 > Y 5 1 


(3*10) 


(3*11) 


After transformation, various terms in the matrix Y^ (eqn. 
3.8) are shown 


4 it «•**,. v 




#*. A. 


Y s = 5 (y 0 +V l + y 2 +y 3 +y 4 +y 5^ 

Y mi = i(y 0 +ay i' a *y2"y3 - “ y 4 +a *y5) 

y m2” y o-“* y l -ay 2 +y 3 - “* y 4~ ay 5 > 

Y m3 = l (y o- y l +Y 2“ y 3 +y 4- y 5 ) 

Y m4 = 5( y 0 -' Iy l- OI *y2 + V3-“ y 4- a * y 5 ) 

Y m5“ |( y o +a * y i- a v 2 -y3-o!*y4+o:y 5 ) 

Now substituting the values of Y°, E° and I' from eqns. (3.8) , 
(3.6) and (3.7) respectively in eqn. (3.2 ) f we get upon simpli- 
fication : 
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For Node N # an additional eqn. may be written by adding 
last six eqn. (3.12) as seen from eqn. (3.4) as follows ; 


C-y 0 ~y 0 -y Q -y 0 -y 0 -y 6y 0 O] 


V, 


p 


V, 


N 


=[I N ] (3.13) 


From eqn. (3.3) we can deduce that 


46 * 

E b “ “V E c 


46 46 46 “ 

•a E , E, = 
a’ d 


*• -* * -x- «■ * 

-E , E = -ocE . E^ = a E 
a e a f a 


Cross multiplying the last six eqns . of (3.12), adding them 
and simplifying the expression gives the following result : 

- y o V l V N " y 0 V2 V N - y o V 3 V N " Wn - Y °Vn - y o V 6 V N 
+6y o V N V N - Y 1 V 1 E 1 ’ “ y l V 2 E a + “ #y l V 3 E a + y l V 4 E t 

+aYl V 5 E* - «* yi V 5 E* + 6 Yl E a E* = S* +S ; + Sc + S*+S* + 4 < 3 - 14 > 


We note from the first seven terms of the eqn. (3.14) that they 
are times eqn. (3.13). If no current is injected through 
the neutral which is the usual condition, 1^=0 and eqn. (3.14) 
simplifies to 


[-y x -aYi **y l Y 1 a Y ± ~ a * Y x 6y£j 


* *)6 * 46 

S +S, +S -hS .+S +S x: 
a b c d e f 

_46 


(3*15) 
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If the neutral is earthed through an admittance y^Q, eqn. 
(3.4) is changed to 


In ■ t 7 +I 8 +I 9 +I 10 +I ir hI 12 +y N0 V N 


(3.16) 


Thus from eqn. (3.13), we can write 


[-y 0 -y 0 -v 0 -y 0 -y 0 -y 0 6y 0 +y NO ] , J “^3 < 3 - 17 > 


Substituting the value of V^I^ in eqn. (3.14), we get with 



I N"° 


[~Yl “ a Yi «*Yi Y± a Yi ~a*y 1 6y 1 ] 


V" 


* * -rc ^ 

S a +S b +S c +S d +S e +S f 


y Npi V N^ 


(3.18) 


Thus the final machine equation (with 1^ = 0) of the six phase 
machine are the first six equations of (3.12), eqn. (3.17) and 
eqns. (3.18) as shown. 


*y x ” a Y 1 a y x y 1 ct Yl -a Yl 6y 1 



In the eqn. (3.19), the term -~g has been neglected as 

E 

is negligible as compared to the “total power in eqn. (3.18). 
Eqn. (3.19) can also be written as 


has been neglected as it 


-*1 0 (6y 0+ y NO ) 1 |v, 


-K- •$£ 

where, S = Sg+S^+S^+S^+Sg+S^ , 


r, * _ * -iT 

[la -a -1 -a a] y x 
[l a -a -1 -a a ] y x 

[min ] 1 y„ 


*■/£■ , or 

S /E a 


(3.2' 
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Eqns. (3.19) or (3.20) is the final eqn. of a six phase machine 
capable of incorporating all the possible unbalances in it. 

This equation is used to form the admittance matrix of six 
phase synchronous generators. 

3.4 MODELLING OF TRANSFORMERS 

The multiphase transformers are normally used to obtain 
6,9,12 and higher phase conversions from three phase systems. 

A six-phase conversion can be obtained from a three-phase 
system by using wye/star, delta/star or some other connections. 
For getting still higher phase conversion, specially constru- 
cted units are used. In this section, symmetrical lattice 
equivalent circuits of two possible transformer connections 
mentioned above and six-phase transformers have been employed 
to obtain the admittance submatrix. The models of other 
schemes may be developed in a similar manner. 

3.4.1 Three Phase/Six Phase Transformers 

The symmetrical lattice equivalent circuit of single 
phase units given in Fig. 2.7 is employed to obtain the model 
of three/six phase transformers. Since most of the time, 
we have got tappings on only one side of the transformer, they 
have been assumed to be t on the primary side and those on 
the secondary side are taken as zero. The parallel windings. of 
three-phase/six phase transformer are assumed to represent an 
equivalent single phase transformer. 
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Consider a wye/star three phase/six phase transformer 
shown in Fig. 3.2. It can be conceptualized to be made of 
three windings [8,9,10] as follows : 

a) a three-phase winding P with terminal a,b,c and 
neutral N 

b) a secondary equivalent three phase winding S having 
terminals 1,£,3 and neutral n 

c) a tertiary equivalent three-phase winding T with terminals 
4,6,2 and neutral n having 180° phase shift with respect 
to the winding S with the neutrals of windings S and T 
joined together to form a common neutral n. 

An equivalent circuit of this transformer can be 
assembled by considering three wye/wye equivalent circuits in 
parallel as shown in Fig. 3.3. If y ps , y pp , Ygy are the 
parameters of the transformer, where y p ,~ in the leakage impe- 
dance referred to the primary side with secondary short cir- 
cuited and tertiary open circuited, and so on, the connection 
table of the transformer can be obtained as shown in Table 
3.1. The three windings P,S and T have turns ratios of aTj_ 

1 and -1 respectively. Since winding T is exactly 180° oppo- 
site in phase to winding S, its turn ratio is negative to that 
of S. Table 3.1 has been derived in a similar manner as done 
earlier in the case of three phase transformers. 


58 


Table 3.1 


Connection table for wye/star three phase/six phase 
transformer with a T = (l+t a ),P T = 1, y t = -1 p.u* 


Admittance 

Between Nodes 

2 

( Yps~^ypT ^ ^ a T 

N-a, N-b, N-c 

YpS +2y ST 

n-1, n-3, n-5 

i 

y P T+2y S T 

n-4, n-6, n-2 

( Vpx^yps ^ 

n-a, n-b, n-c ' 

W“t 

a-1, b-3, c-5 

-y ps /«T ; 

N-l, N-3, N-5 

~Yp T / a T ' 

a -4, b-6, c-2 

y pT /a T 

N-4, N-6, N-2 

- y ST 

1-4, 3-6, 5-2 


From Table 3*1, the nodal admittance matrix of the above trans- 
former can be obtained as shown : 
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abcl2345 6 



y 3 ~ (YpT +Y si) * Y 4 " Yps/ a x 

y 5 = yp T / a T ; y 6 = yst 

The nodal admittance matrix of a delta/star transformer shown in 
Fig* 3*4 can also be obtained in a similar-manner. In this 
case, also, we can visualize the transformer as consisting of 
, three windings with the change from the case of wye/star trans- 
former that primary is a delta connected three phase winding 
Mth turns ratio a T = f3 (l+t a ), where t a is the off-nominal 
pping on the primary side. Thus, the equivalent circuit 




FIG. 3.2 A THRE PHASE/SIX PHASE , WYE/STAR TRANSFORMER 
OF TURNS RATIO <V T : 1 



FIG 3 3 REPRESENTATION OF A WYE/ STAR TRANSFORMER FOR 
OBTAINING AN EQUIVALENT CIRCUIT. 


61 


of three phase/six phase delta/star transformer is obtained by 
paralleling two delta/wye and one wye/wye equivalent circuits 
in turn as shown in Fig. 3.5. Table 3.2 gives the connections 
showing admittances between various nodes with the same assump- 
tion that parallel windings represent single phase units. 

Table 3.2 

Connection table for delta/star three/six phase 
transformer with a T = f3(l+t a ), p =0L, y T = -1 p.u. 


Admittance 

Between Nodes 1 

( yps+Yp T )/aj 

a-b, b-c, c-a 1 

» ; 

( Yps +2 YST^ 

i 

1-n, 3-n, 5-n ! 

( Ypx’^YsT^ 

4-n, 6-n, 2-n I 

y ps /a T 

1-c, 3-a, 5-b 

-Y ps / a T 

1-b, 3-c, 5-a 

'~ y PT^ a T 

4-c, 6-a, 2-b 

y p i/a T J 

4-b, 6~c, 2-a 

” y ST 

i 

1-4, 3-6, 5-2 
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The nodal admittance matrix obtained from Table 3.2 is given 
below ; 


a b c 123456 


a 

b 

c 

Y 1 Y 2 y 2 

Y 2 Y 1 Y 2 

Y 2 Y 2 Y 1 

y 3 y 4 -^4 -y 3 

_y 4 - y 3 y 4 y 3 

y 4 -y, ; -y 3 y 3 . i 

Y = 

TR2 


(3.22): 

1 

-y 4 y 4 

~ Y 5 

2 

y 3 -*3 

“ Y 5 

3 

y 4 -y 4 

" y 5 ' 

4 

y 4 -y 3 

; 

“ Y 5 

5 

~ Y 4 Y 3 

“ y 5 

6 

~ Y 3 y 3 

-y 5 


a r) 

where, y x = 2(y ps +y pT )/a T , y 2 = (y ps +y pT )/a T , 

y 3 ~ y PT^ a T' y 4 = Yps/ a T» y 5 “ y ST 

3.4.2 Six Phase Transformers 

Consider a star/star six phase transformer shown in Fig. 
3.6. Assuming that parallel windings represent single units, 




FIG. 3 4 A THREE-PHASE/SIX Pl'ASE DELTA/STAR TRANSFORMER OF 
TURNS RATIO <X T :1 



FIG. 3.5 REPRESENTATION OF A' DELTA/STAR TRANSFORMER FOR 
OBTAINING AN EQUIVALENT CIRCUIT . 
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admittances between various nodes are obtained as shown in 
Table 3.3. 


Table 3.3 

Connection table for a six phase star/star -transformer 
with a T = (l+t a ), p T « 1 p, u . 


I 1 "" ■ 1 111 "" n " ' ■ - Ji 

Admittance 

Between Nodes 

y/«x 

N-a, N-b, N-c, N-d, N-e, N-f 

y 

n-1, n-2, n-3, n-4, n-5, n-6 

y/a T 

a-1, b-2, c-3, d-4, e-5, f-6 

-y/«x 

n-a, n-b, n-c, n-d, n-e, n-f 

-y/a T 

N-l, N-2, N-3, N-4, N-5, N-6 

6y/ Oj 

N-n 


The nodal admittance matrix for such transformer as obtained 
from Table 3.3 is 


a 


f 1 . • • • 6 


TR3 


f 

1 

♦ 

♦ 

6 


Y T1 

Y T2 

_ 

‘ ‘ 

Y T2 

Y T3 

< _ 



(3.23) 
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where , 

y t1 = <*iag. Cl 1 1 1 1 l](y/c^) 

Y t2 = diag. [l 1 1 1 1 l](-y/a T ) 

Yj 3 = diag. [l 1 1 1 1 l]y 

where y is the series leakage admittance in p.u. of the 

transformer. 

In deriving the representations for transformers, the 
effect of magnetising core and saturation are neglected. It 
has been found out that the models derived above are quite 
precise for unbalanced studies of the system. The representa- 
tions for other types of transformer can also be obtained in 
the manner illustrated above. 

3.5 MODELLING OF TRANSMISSION LINES 

It is necessary to derive the models of transmission 
line which forms an integral part of the power systems. In 
this section, the representation of multiphase lines in terms 
of phase Impedance matrix, Tt-circuit representation and ABCD 
parameters have been developed. Three phase equivalent repre- 
sentation of a six-phase line and vice versa have also been 
done. 
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a) Phase Impedance Matrix : 

A six phase short ^transposed transmission line processing 
cyclic symmetry is described by its phase impedance matrix as 
shown 


24) 


This matrix can be fully diagonalized by using symmetrical 
component six phase transformations. However, medium and long 
lines cannot be representation by (3.24). 

b) it -Circuit Representation : 

The three phase nominal and equivalent n-circuit represen 
tation for three phase systems can be extended to six phase 
systems as well to adequately represent medium and long lines 

as shown : 


I 6 

x pS 

S3 

Y^+ *r Y^v, 

*p + 2 1 sh 

-Y 6 

P 


J 

(3.25) 

x pR 

* 

y6 

P 

■ 



» \ 


4 


Z 11 

Z 12 

Z 13 

Z 14 

Z 15 

Z 16 

hi 

h.2 

h.3 

h.4 

Z 25 

h>6 

Z 31 

Z 32 

Z 33 

h>4 

Z 35 

Z 36 

Z 41 

Z 42 

Z 43 

Z 44 

Z 45 

Z 46 

hi 

Z 52 

Z 53 

Z 54 

Z 55 

Z 56 

Z 61 

Z 62 

Z 63 

Z 64 

Z 65 

Z 66 


(3, 
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The derivation of eqn. (3.25) can bo done on the same lines 

as for three phase transmission lines. In eqn. (3.25), the 

subscripts S and R denote sending and receiving ends respecti- 
£ 

vely, Y° is a (6x6) phase admittance matrix obtained from 
H 

the relation (3.9). 


c) ABCD-Parameters Representation : 

The representation of six phase transmission lines using 
ABCD parameters can be written as follows : 






I 


6 

pS 
**** ' 




(3.26) 


d) Equivalent Three Phase Representation of a Six-Phase Line : 

If we are interested in carrying out the analysis of a 
composite three phase six phase system, the interest of investi- 
gation being on the three phase part of it, we have to formulate 
an equivalent three phase representation of the six phase part 
of the system. For example, consider a system having a six 
phase transmission line connected between three phase buses S 
and R via three phase/six phase transformers TR X and TR^, as 
shown in Fig. 3.7. The equation of the transmission line 
connected between buses S’ and R' in impedance form is written 


as 


- v t- ■ 


(3.27) 
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Assuming the transformers TR^ and to be wye/star, 
three phase/six phase transformers, we can write the expre- 
ssions for three phase voltages Vg and v2 in terms of six 
phase voltages V^,, V^, and current 1^ [10] as follows ; 

v| = 5 ^ + k nT I p (3.28) 

V R = 2 nTv r- - 2 nT t-hl 4 (3,29) 


where [ Z-^] and are (6x6) diagonal matrices representing 

leakage impedances of the transformers TR^ and TR 2 and [ll] 


N = 


10 0 
0 0-1 
010 
-10 0 


0 0 1 
0-10 


In addition, the terminal relationship for a wye/star transfor' 
mer shown in Fig* 3.2 may be written as 


v£ = vjj, - [Z] Ip 

p p F 

Subtracting eqn. (3.29) from eqn. (3.28), we obtain 


(3.30) 


v2 - vj* * k n t (v£, - 4* } + 5 nT t z i + h) 1 


6 


R* 


(3.31) 
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Using eqn. (3.27) and the relationship I 6 

P 

v i - ^ - 2 NT [Zp + + %] »i p 3 

Hence, from eqn. (3.32), we obtain 

^p,eq ~ 2 ^ + ^1 + ^2^ (3.33) 

The equivalent three phase admittance matrix can either be 

obtained by inverting zj or derived as follows : 

p , eq 

From eqn. (3.27) we get 

Ip = Vp [v|, - V®,] (3.34) 

Substituting the values of V~, and V~, from eqns. (3.28) and 
(3.29) respectively in eqn. (3.34), we obtain 

Ip = Y p N[vf - v|] + [Z 1+Z2 ] Nip (3.35) 

We also have [ll] 

if = £ N T i| (3.36) 

Using eqn. (3.35) we get 

Y p,eq = ^ U+ i N ^ 1 $ 0.37] 

where U is a (3x3) identity matrix • 


= Nip [11], we get 
(3.32) 
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In a similar fashion, we can derive the three phase 
equivalent representations when other types of transformers 
are used. 


The three phase equivalent representation based upon 
A3CD parameters can also be derived. From Fig. 3.7 we get 
the terminal relationships at bus bars S, S' , R and R‘ as 
shown 


V 


pS 
3 

:p s J 


5 nT 


0 


k NT £ z lJ 


1 m t 

N 


V' 


pS' 

[ 6 

L pS' 


(3.38) 


and 


V 6 - 
V pR’ 


'pR’ 


N 


[Z^N 


N 


v: 


pR 

[3 

pR 


Using eqns. (3.26), (3.38) and (3.39) we can write 


(3.39) 




■pS 




2 N<r 


0 


k N T [Z x ] 


i-N T 


5 


B 


D 


N [Z^N 


0 N 




P R J 
(3.40) 



. * 

-I 3 1 
A pSJ 

' & 1 

'* 


A* 

^1 eq 


i 3 1 
1 eq 


4 


eq 


d; 


1 eq 


V 3 ’ 

V pR 


PR 


(3.41) 



_ - - -J 2 £ 

Six phase transformei] 3 


f Rj 


line 


R / 


tr 2 


R 


FIG. 3.7 A SIX -PHASE TRANSFORMER LINE CONNECTED BETWEEN 
TWO THREE PHASE BUSES VIA 3-0/ 8-0 TRANSFORMERS 
TR t AND TRj 4 , 
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FIG. 3.8 A THREE PHASE TRANSMISSION LINE CONNECTED BETWEEN 
TWO 6-0 BUSES THROUGH 6-0 /3-0 TRANSFORMERS TRi 8, 
TR 2 . 

*3 j fe Jp ^ lh 

3-0 3-0Z6-0 6-0 

But Transformer Load 

PIG. 3.9 A SIX PHASE LOAD CONNECTED TO THREE PHASE BUSBAR 


I t «Lj |— i 

I 9 r 
6-0 6-0 /3-0 

Bus Transformer 


FIG. 3.10 A THREE PHASE LOAD CONNECTED TO SIX PHASE BUSBAR 



3-0 

Load 
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Noting that A 6 = D 5 and A. B, C, D , [zj and [ 22 ] are diagonal 
matrices, eqn, (3.41) after simplification gives 

A?eq £ N T A 6 N + £ N T [z 1 ] C 6 N 

D?eq - £ N T B^N + £ + £ N T [Z 1 ]C 6 [Z 2 ]N 

C?eq « N^C^N 

D?eq * ^ N T D 6 N + £ N T C 6 [Z;^]!'! 

e) Equivalent six phase representation of a three phase line : 

If we are interested in the analysis of the six phase part 
of a composite three phase/six phase power system, an equivalent 
six phase representation of the three phase part of it becomes 
essential. However, the analysis on equivalent six phase basis 
may become tedious because of prohibitive memory requirements. 

Consider a system having a three phase line connected to 
six phase buses via six phase/three phase transformers TR^ and 
TR 2 as shown in Fig, 3.8. The six phase equivalent impedance 
and admittance may be derived [ll] in a similar fashion as done 
earlier in subhead (d). The final results are : 

Z p f @q * | Zj* N T + (3.42) 

Y p,eq * I « Y p ^[Zi+Z^r 1 5 N N T 


(3.43) 


73 


Noting that = D ^ and A, B, C, D , [Z^] and [z 2 ] are diagonal 
matrices, eqn. (3.41) after simplification gives 

A?eq = j N T A 6 N + ^ N T [Z 1 ] C 6 N 

B?eq = ^ N T B 6 N + ^ N T A 6 [Z 1 +Z 2 ]N + ^ N T [Z i ]C 6 [Z 2 ]N 
C?eq = *r N T C 6 N 

D?eq = ^ N T D 6 N + ^ N T C 6 [Z^H 


e) Equivalent six phase representation of a three phase line : 

If we are interested in the analysis of the six phase part 
of a composite three phase/six phase power system, an equivalent 
six phase representation of the three phase part of it becomes 
essential. However, the analysis on equivalent six phase basis 
may become tedious because of prohibitive memory requirements. 

Consider a system having a three phase line connected to 
six phase buses via six phase/three phase transformers TR^ and 
TR 2 as shown in Fig. 3.3. The six phase equivalent impedance 
and admittance may be derived [ll] in a similar fashion as done 
earlier in subhead (d). The final results are : 

Z p,eq = 5 N 4 N T + t z l +z d (3>42) , 

Yp >eq = [U+ 2 N Yp N T [Zi + z 2 ]]" 1 5 N Y 3 N T (3.43) 
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where U is a (6x6) order identity matrix. 

Similarly, the equivalent six phase ABCD parameters of the 
transmission line may be derived as under : 

A?eq = 5 N A V + § [ Zl ]NON T 

B?eq = j NB 3 N T + ^ NA 3 N T [Z 2 1+ \ [Z 1 ]ND 3 N T + 

+ \ [z 1 ]nc 3 n t [z 2 3 

C^eq = j NC 3 N T 

D^eq = | ND 3 N T + ^ NcV^] 

3.6 MODELLING OF SIX PHASE LOADS 

Multiphase loads may be modelled as a constant impedance/ 
admittance to ground. Here again, we have to obtain three phase 
equivalent representation of a six phase load or vice versa 
depending upon whether our interest of analysis lies on the 
three phase part or six phase part of the composite system. 

Consider a six phase load connected via a three phase/ 
six phase wye/star transformer to a three phase bus as shown 
in Fig. 3.9. 

For the six phase load, we get 

s= 1 ^ 
v p no p 


(3.44) 
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We also get voltage V 3 as 

P 


V 3 

p 

= \ N T v£ + | n t [z] ip 

(3.45) 

or 



V 3 

p 

= \ n t z l 6 ni 3 + \ n t [z]ni 3 

(3.46) 

Therefore, 



z 3 

L,eq 

= 5 n t [z* + z] n 

(3.47) 


Similarly, we can derive the equivalent three phase load 
admi ttance matrix as shown below : 


Y L,eq ~ + 2 N T Y^[z]N ]’” 1 N T y£n] (3.48) 

Now, consider a three phase load connected to a six phase bus 
through a six phase/three phase star/wye transformer as shown 
in Fig. 3.10. The load can be represented as 


v: 


= z 3 i 3 

z L i p 


(3.49) 


We obtain voltage V” as 


v p = 5 N Zl nT]1 p + ^ 1 


(3.50) 


or 


Z L,eq =5 WZ l' nT+Z 


(3.51) 

Similarly, we can derive the six phase equivalent admittance 
in phase form as shown : 

yf = [U+ i NY^Z]” 1 [i NY^N T 3 (3.52) 
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3.7 SAMPLE NETWORK 

The mixed three phase and six phase load flow analysis of 
the sample network shown in Fig. 3.11 was carried out. The 
sample system is the same as that of Fig. 2.14 except that the 
three double circuit three phase lines are replaced by six 
phase lines. This can be achieved by installing two three- 
phase/six-phase wye-star transformers at the sending end and 
at the receiving end respectively of the six phase line. The 
bus numbering sequence with such an arrangement is shown 
in Fig. 3.11. 

The effect of converting the double circuit three phase 
line to a six phase line has been studied. For the purpose of 
study, each double circuit three phase line was converted to a 
six phase line one by one. Then, two double circuit three phase 
lines were converted to six phase lines simultaneously and 
finally, all the three double circuit three phase lines were 
converted to six phase lines. The necessary data are given in 
Table 3.3. Gauss-Seidel iterative technique was used to 
conduct the load flow analysis. The algorithm used was the same 
as that used in the three phase load flow except that the 
admittance submatrices for the six phase lines and three phase/ 
six phase transformers were formed separately and the acquittance 
matrix of the system modified accordingly. The flow chart for 
developing these submatrices is given in the Appendix. 
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In the first instance, the load flow analysis of the 
balanced system were conducted for all the cases mentioned above 
After that, the system was made unbalanced with the introduction 
of the same unbalances which were simulated in the case of three 
phase load flow (Chapter II) in order to study their effect on 
the system when it contains six phase lines. The results 
obtained from the balanced and unbalanced load flow studies with 
three six phase lines in the system are shown in Tables 3 . 4 — 3 .5. 

3.8 CONCLUSIONS 

The relative merits and demerits of multiphase systems over 
the conventional three phase systems have been discussed. The 
detailed mathematical modelling of the six phase elements such 
as machines, transformers and transmission line has been deve- 
loped. Three phase equivalent representation of a six phase 
line and vice versa have also been developed. The following 
conclusions are drawn from the mixed three phase six phase load 
flow analysis conducted and results given in Tables 3.4— 3*5. 

a) When only one double circuit three phase line is converted 
to a six phase line, the voltage magnitudes at the two buses 
connecting the six phase line are improved while the voltages 
at the other buses remain nearly unchanged. 

b) The improvement in the voltage magnitudes is more when two 
double circuit three phase lines are converted to six phase 
lines and the results are best when all the three double circuit 
three phase lines are six phase lines. 
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c) The voltage angles at the buses near the generators are 
reduced slightly. 

d) As a consequence of reduction of voltage angles, the 
reactive power flows from the generators are reduced slightly. 

e) The real power flow from the slack generator increases 
slightly . 

f) It was found, that, with the same unbalances introduced 
as in the three phase case, the degree of unbalance in the 
case of mixed three phase six phase load flow remains same. 

In fact, there is a slight improvement in the voltage magni- 
tudes unbalancing. 

The above para give the advantages of the six phase 
transmission system over the conventional three phase systems. 



T2 U-6) (31-33) 

* 
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P* 180 MW 
IVI = 1.04 V 


FIG, 3.11 THE SAMPLE NETWORK WITH BUS NUMBERING SEQUENCE 
FOR MIXED 3-0 AND 5-0 LOAD FLOW STUDY. 
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Table 3.3 

Circuit characteristics of the network of Fig. 3,11 on 
100 MVA/phase base and rate KV 


~rom To 

Bus Bus 

12 3 

... i J i 

4,5 sequence 


Zero sequence 



R 

X 

B 

R o 

X 

o 

B o 

L) 6-0 Lines : 







(40-45) 

(46-51) 

0.0125 

0.030 

0.01875 

0.03875 

0.145 

0.01087' 

(52-57) 

(58-63) 

0.00625 

0.025 

0.00938 

0.01934 

0.0725 

0.0054 

(64-69) 

(70-75) 

0.0083 

0.033 

0.01245 

0.0258 

0.0966 

0.0072 


-) 3-0/6-0 Transformers : 

The transformers T^-Tg are wye-star solidly earthed with 
YpS = Ypy = Yst = *"Jl2.5 

(Rest of the data is same as that in Table 2.4) 


81 


Table 3.4 

Voltage solutions for mixed three phase six phase load 
flow for the network of Fig. 3.11 


Bus 

Balanced Case 

Unbalanced Case 

No. 

I vl 

6 

|v| 

6 


in p*u. 

in deg. 

in p.u. 

in deg. 

1 

1.043 

-2.731 

1.039 

-5.111 

2 

1.043 

-120.799 

1.034 

-124.712 

3 

1.043 

11S.335 

1.044 

115.446 

4 

1.027 

-1.449 

1.032 

89.588 

5 

1.028 

*1X9.504 

1.016 

-28.871 

6 

1.028 

+119.575 

1.040 

-148.857 

7 

1.006 

-3.497 

1.003 

-5.995 

8 

1.007 

-121.554 

0.998 

-125.479 

9 

1.006 

117.525 

1.009 

114.616 

10 

1.031 

82.518 

1.024 

82.862 

11 

1.032 

-37.411 

1.034 

-37.700 

12 

1.031 

-157.466 

1.017 

-157.840 

13 

1.026 

81.927 

1.020 

82.257 

14 

1.027 

-38.002 

1.032 

-38.309 

15 

1.027 

-158.056 

1.014 

-158.412 

16 

1.024 

83.599 

1.016 

84.913 

17 

1.024 

-36.332 

1.020 

-36.902 

18 

1.025 

-156.365 

0.988 

-157.115 

19 

0.993 

80.377 

0.987 

81.016 

20 

0.993 

-39.558 

1.000 

-39.921 

21 

0.994 

-159.588 

0.973 

-160.294 

22 

0.997 

80.821 

0.991 

81.382 

23 

0.997 

-39.109 

1.003 

-39.458 

24 

0.997 

-159.145 

0.979 

-159.789 

25 

1.022 

81.784 

1.017 

82.182 

26 

1.023 

-38.145 

1.028 

-38.454 

27 

1.023 

-158.198 

1.009 

-158.631 

28 

1.070 

0.000 

1.070 

0.000 

29 

1.070 

-120.072 

1.070 

-120.072 

30 

1.070 

120.072 

1.070 

120.072 

31 

1.040 

2.995 

1.040 

96.102 

32 

1.041 

-117.062 

1.041 

-23.955 ' 

33 

1.041 

123.003 

1.041 

-143.968 i 

34 

1.020 

-1.729 

1.020 

-1.829 , 

35 

1.021 

-121.787 

1.021 

-121.886 j 

36 

1.021 

118.278 

1.021 

118.179 j 

37 

0.000 

0.000 

0.000 

0.000 ! 

38 

0.000 

0.000 

0.000 

0.000 

39 

0.000 

0.000 

0.000 

0.000 

40 

1.782 

82.432 

1.773 

82.774 ! 
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Bus No. 

Balanced Case ’ 

Unbalanced Case 


l vi 

0 

M 

9 


in p.u* 

in deg. 

in p*u. 

in deg. 

41 

1.782 

22.449 

1.761 

22.076 

42 

1.782 

-37.496 

1.791 

-37.787 

43 

1.782 

-97.567 

1.773 

-97.226 

44 

1.783 

-157.551 

1.761 

-157.923 

: 45 

1.782 

142.504 

1.791 

142.214 

46 

1.778 

82.012 

1.768 

82.346 

47 

1.778 

22.029 

1.756 

21.672 

48 

1.778 

-37.917 

1.788 

-38.222 

49 

1.778 

-97.988 

1.768 

-97.654 

50 

1.778 

-157.970 

1.756 

-158.327 

51 

1.778 

142.084 

1.788 

141.779 

52 

1.726 

80.734 

1.716 

81.308 

53 

1.726 

20.768 

1.694 

20.109 

54 

1.726 

-39.196 

1.737 

-39.548 

55 

1.726 

-99.265 

1.716 

-98.691 

56 

1.726 

-159.231 

1.694 

-159.890 

57 

1.726 

140.805 

1.737 

140.453 

58 

1.721 

80.464 

1.711 

81.089 

59 

1.722 

20.499 

1.687 

19.808 

60 

1.721 

-39.466 

1.733 

-39.831 

' 61 

1.721 

-99.535 

1.711 

-98.910 

: 62 

1.722 

-159.501 

1.687 

-160.191 

: 63 

1.721 

140.535 

1.733 

140.170 

; 64 

1.780 

82.348 

1.771 

82.703 

1 65 

1.781 

22.364 

1.759 

21.975 

: 66 

1.781 

-37.581 

1.789 

-37.874 

67 

1.780 

-97.652 

1.771 

-97.297 

68 

1.781 

-157.636 

1.759 

-158.025 

69 

1.781 

142.419 

1.789 

142.127 

70 

1.772 

81.955 

1.763 

82.342 

71 

1.773 

21.973 

1.750 

21.556 

72 

1.773 

-37.974 

1.782 

-38.279 

73 

1.772 

-98.045 

1.763 

- 97.657 

74 

1.773 

-158.027 

1.750 

-158.444 

75 

1.773 

142.027 

1.782 

141.721 






83 


Table 3.5 
Power Flows 


To 

Bus 

From 

Bus 

Balanced Case 

Unbalanced Case 

Real Power 
(MW) 

Reactive 

(MVAR) 

Real Power 
(MW) 

Reactive Power 

(mvar) 

1 

28 

95.102 

31.551 

102.490 

38.926 

2 

29 

95.032 

31.350 

92.591 

43.812 

3 

30 

95.229 

31.552 

93.213 

32.891 

' 4 

31 

59.675 

10.536 

70.863 

9.017 

5 

32 

59.759 

10.715 

53.368 

17.509 

6 

38 

59.546 

10.642 

54.743 

-1.960 

7 

34 

39.687 

9.802 

43.763 

11.793 

8 

35 

39.749 

9.968 

37.576 

15.258 

9 

36 

39.564' 

9.899 

37.663 

8.124 

25 

22 

-11.050 

-11.561 

-9.725 

-12.321 

’26 

23 

-11.100 

-11.716 

-11.790 

-11.070 

27 

24 

-10.948 

-11.686 

-12.440 

-12.985 

25 

13 

0.846 

-2.715 

0.609 

-2.611 

26 

14 

0.849 

-2.710 

0.894 

-2.864 

27 

15 

0.841 

-2.712 

1.081 

-2,496 

16 

13 

- 4.727 

-3.718 

-7.290 

-2.103 

17 

14 

-4.707 

-3.656 

-3.968 

-4.851 

: 18 

15 

- 4.775 

-3.671 

-2.791 

- 0.273 

25 

16 

7.965 

-4.792 

11.650 

-6.681 

; 26 

17 

7.937 

-4.883 

6.869 

-3.424 

■ 27 

18 

8.034 

-4.860 

5.276 

-9.352 

19 

16 

30.865 

7.291 

35.958 

4.766 

20 

17 

30.854 

7.271 

29.812 

9.287 

21 

18 

30.892 

7.283 

27.061 

-0.026 

■ 40 

46 

-13.979 

-1.612 

-14.331 

-2.599 

41 

47 

-13.967 

-1.637 

-13.423 

-2.740 

42 

48 

-13.989 

-1.637 

-14.442 

-0.928 

43 

49 

-13.979 

-1.611 

-14.330 

-2.599 

44 

50 

-13.967 

-1.636 

-13.423 

-2.739 

45 

51 

-13.988 

-1.636 

-14.441 

-0.927 

64 

70 

- 27.879 

-9.160 

-25.827 

-10.356 

65 

71 

-27.818 

-9.247 

-29.648 

-12.017 
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To 

Prom 

Balanced 

Case 

Unbalanc ed 

Case 

■Duo 

JDU.S 

Real Power 
(MW) 

Reactive 

Power 

(MVAR) 

Real Power 
(MW) 

Reactive : 

Power i 

(MVAR) 

66 

72 

-27.910 

-9.267 

-28.772 

-8.187 

67 

73 

-27.879 

-9.159 

-25.826 

-10.356 

68 

74 

-27.817 

-9.247 

-29.648 

-12.016 

69 

75 

-27.910 

-9.267 

-28.771 

-8.186 

52 

58 

-13.347 

-4.118 

-11.197 

-5.378 

53 

59 

-13.345 

-4.127 

-15.037 

-7.184 

54 

60 

-13.349 

-4.123 

-13.919 

-3.095 

. 55 

61 

-13.347 

-4.118 

-11.197 

-5.378 

56 

62 

-13.345 

-4.126 

-15.037 

-7.184 

57 

63 

-13.349 

-4.123 

-13.919 

-3.095 



CHAPTER 4 


TWELVE PHASE POWER SYSTEM: MODELLING 

4.1 INTRODUCTION 

As the demand for electric power is increasing at an 

exponential rate, the emphasis of the scientists $nd power 

* 

* • - 

system engineers has been to develop feasible dlternatives of 
transmitting more power from generating stations to the load 
centres. One of the striking alternatives, as pointed out 
earlier, is the use of high phase order transmission. 

A considerable amount of work has been done on the 
feasibility of six phase systems, yet the topic of twelve 
phase systems seems to have been left untouched, in this 
respect. Although the idea of using twelve phase systems is 
in its infancy, one should not wonder if it becomes a 
reality in future. To make twelve phase systems operative, 
careful planning of such system is required. Planning studies 
can be done only if suitable mathematical modelling of the 
twelve phase systems is available. 

This chapter presents the mathematical modelling of the 
different elements present in a twelve-phase power system 
network. In the first part, the twelve phase synchronous 
machine equation incorporating all the unbalances has been 
developed. After that, the nodal relationships between 
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current and voltage for a threerphase/twelve phase transformer 
has been derived. This is followed by the presentation of 
admittance matrix of a twelve phase transmission line based 
upon the phase impedance matrix, it circuit representation and 
ABCD parameters. Three phase equivalent representation of a 
twelve phase line and vice versa have also been derived. Fina- 
lly, equations have been developed to model twelve-phase loads* 

4.2 MODELLING OF TWELVE -PHASE GENERATOR 

Although the feasibility of a twelve-phase generator is 
yet to be established (even six-phase generators have not yet 
come into existance) , a model similar to that of three and 
six-phase synchronous generators has been developed in this 
section. 


The schematic representation of a general twelve-phase 
element is given in Fig. 4.1. The figure represents a star- 
connected twelve-phase synchronous machine if E a ,E b , ...» 
are a balanced set of internal voltages and nodes 13, 14,..., 24 
are joined together to form a star point. Hence, a twelve 
phase synchronous generator is represented by a balanced set 
of internal voltages behind proper reactances. The equation 
of the element on Fig. 4.1 in the impedance form is 


^ n 12 _ -12 12 12 
V-, + E- = Z “I n , + V n 


r 12 

P' 


*P 


p x p 1 


( 4 . 1 ) 




FIG. 4.1 SCHEMATIC REPRESENTATION OF A GENERAL 12-0 ELEMENT. 


© © 
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where 


V 12 = 
p* 

^ V 13 

V 14 *•* 

v 24 ] t 

p 12 _ 
E P ~ 

[E a 

E b 

e x ] t 

V 12 ' = 
P 

[Vi 

v 2 ... 

v i2 ] t 

T 12 
v - 

^13 

I 14 **• 


2= 

t 1 . 

I b • •• 



12 

and z p » (12x12) impedance matrix of the 12-phase element 

* [v^r 1 


From eqn. (4,1), we can write the nodal equations of Fig. 4.1 
for the voltage and current directions shown as follows : 


I 


12 

P 


y 12 12 
P Y p 



where, I* 2 = [x 1 I 2 ... 


_yl2 

P 


V 2 1 

p 

CM 
r-i i 
>* 


V 12 

p* 


E 12 



p 



L— 

Xi,] 1 

" 



(4.2) 


For a balanced set of internal voltages, 

— 12 r , * 2 * 2 * O-H- O T 

E p * L 1 a a -a -or -a -1 -a -a 2 j a 2 a ] T E a (4.3) 

u 

where, a = e j 2 */l2 = o.866+j0.5 * 12 fl = i£jo° 
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If the nodes 13-24 are joined together, then 


I N = X 13 +I 14 + + X 24 


V N ~ V 13 ~ V 14 = ’ * = V 24 


(4*4) 

(4.5) 


Thus internal voltages are given by, 


E p “ ^ E a +V N' E b +V N» * **» e 1 +v N^ 


(4.6) 


and the current by 


* * 

_3_- . r _ s b . 

< W* * b " ( V E b )* ’ 


(VHP 


(4.7) 


where S a ,S^, ..., are complex powers for phases a,b,...,l 
respectively. 

The (12x12) admittance matrix in eqn. (4.2) is given 


s 

Y ml 

Y m2 

Y m3 '** 

Y mll 

mil 

Y S 

Y ml 

Y m2 •** 

Y mlO 

mlO 

Y mll 

Y , 

• 

Y ml * ' * 

Y m9 

ml 

Y m2 

« 

• 

Y m3 

Y m4 **• 

Y s 


(4.8) 
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,12 . 


The expression for Yp is obtained after using symmetrical coropo- 


1 o 

nent twelve-phase transformation matrix T [24] 

5 


i *e 


• * 


y 12 = [t 12 ][y 12 ][T 12 ] T 

p 1 s JL comp JLA s J 


(4.9) 


where. 


t 12 = 1 _ 
s -yi2 


and Y 


12 

comp 


1 

1 

i i 

1 

1 

i 

1 

1 

1 

1 

1 

1 

a 

a 2 j 

2* * 

-a -a 

-i 

-a 

-a 2 


a 2 * 

* 

a 

1 

a 2 . 

„2* . 
-a -l 

2 
— oc 

„2* 

a 

i 

a 2 

-a 2 * 

-1 

-a 2 

„2* 
c t 

1 

j ■ 

-1 -j 

1 

3 

-i 

* 

-J 

1 

j 

-1 

-j 

1 

2* 2 , 
-a -a 1 

2* 2 
-a -a 

i 

2* 

-ar 

-a 2 

1 

-a 2 * 

-a 2 

1 

* 

--a 

2* . 
a J 

-a 2 

a 

-a 

* 

a 

-a 2 * 


a 2 

-<x 

1 

-1 

1 -1 

1 

1 

-1 

1 

-1 

-1 

1 


1 

-a 

J1 . 
a -j 

-a 2 * 

* 

a 

-1 

a 

„2 

—a 

♦ 

3 

-a 2 * 

-a 

1 

-a 2 • 

„2* . 
-a 1 

-a 2 

2* 

-a 

1 

-a 2 

2* 
— oc 

1 

-a 2 

2* 

-a 

1 

-j 

-1 j 

1 

-j 

-1 

J 

1 


-1 

j 

1 

2* 

a 

J2 , 

-a -l 

2* 

-a 

a 2 

1 

2* 

a 

^2 
— ct 

-1 

-a 2 * 

a 2 

1 

* 

a 

2* 

a -j 

2. 

— cc 

-a 

-1 

* 

-a 

2* 

-a 

j 

a 2 

a 

- 

=s 

diag* 

[y 0 . 

Yl, Y 2 

» * * 

1 • t 

Yll ] 






(4.11) 


y 0 ,Yi,y2 »*»« are 'twelve-phase sequence admittances. 

After using eqn. (4.9), various elements of the eqn. (4.8) are 
obtained as given below : 
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Y s ( y o + yi + y2 + y3 + y4 + y5+y6 + y7+y8 + y9+y 10 + yii) 

Y ml = 12 (y 0 +a yi +a2 y 2 + 3 y 3 -^y 4 -o: # y 5 -y 6 -ay 7 -a 2 yg-jy 9 +a^y 1 o+^y 11 ] 
Y m2 t? (y 0 +a vi-« Y 2 ~y 3 ~ a y 4 +a y5 + y6 +a Y7“ a2 y 8 

2 o* 

-y 9 ~a y 10+ a y u ) 

Y m3 = 12 < y 0+ j yi-y 2 -o y 3 +y 4 +3 y 5 -y 6 - j y 7 +y 8 + j y 9 -y 10 ~ j y u ) 

Y m4 = 12 (y 0 -a 2 yi“ a2 y2 + y 3 ’- a2 ''y4- a2 y5 + y6'- a2 ''y7“ K2 y8 + Y9^ 2 *yio 

- a 2 y n ) 

Y ms~ 12 ty 0 ~ a yi +a y 2 + Jy 3 - a y4 +a Y5-y6 +a y 7 ~ a y 8 -jy 9 

+ “ 2y io _0Iy n) 


Y m6~ i2 ^yo“yi + y2-y3 +y 4-y5 + y6’-y7 + y8“y9 + yio-*yn) 

v — 1 / . 2 . 2 *** o 

Y m7~ 12 (y<>~ a Yl +a y 2 -jy 3 -a y 4 +a y 5 -a y 6 +ay ? -a y 8 +jy 9 

2 * * 

+ a y 10 -o y u ) 

Y m 8 = T 5 ( y 0 - a 2y l- a2 * y 2 +y 3 -c ‘ 2y 4' a2,ty 5 +y 6 - “ 2y 7- ct2 * y 8 

2 2* v 

+ y 9 -a y 10 -o y u ) 
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Y m9 = 12 (y o -^l->'2 + jY3 + y4-^5-y6 + ^7 + V8-^9' Y 10 + ^ll ) 

1 / 2* 2 2* 2 2* 2 
Y mlO = T2 (y o +ct y l"“ y 2 -y 3- a y 4 +a y 5 +y 6 + “ y 7 _a y 8 

-yg-a^Yio+^yu) 

, I . * 2* 2 * 2* 

Y mll = 12 ^ V a yi +a y2~jy3~ a ^4“ a y5^6“ a *7“® ^8 

+3Y 9 +cc 2 Y 10 +ay 1 ) 

Substituting the values of E^, 1^? and from eqns. (4,6), 
(4.7) and (4.8) respectively in eqn. (4.2), we obtain after 
simplification and collection of terms 
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y l 



a ■Y^ 

a z y 1 



JYl 
— a 2 *y. 


V N 

1 

% 

•c 


Tl 

1 

C 

-“Yi 


E» 

a 

-a2 yi 



-3 y i 

a^Yi 



a y l 


_ - 


s b/<V E b>* 


s*/(v n+ Ei)*J 

(4.12) 


By adding the last twelve eqns. of (4.12), we get an additional 
equation for node N as 


= 


C“Y 0 ~Y 0 ~Y 0 -y 0 -Y 0 ~Y 0 -y 0 -Y 0 -y 0 -y 0 -y Q y Q 12y 0 0] 

Mm -*• 

r 12 


v: 


v, 


N 


(4.13) 


It can be deduced by taking the conjugate of eqn. (4*3) that 


1 2 r - 2.2**. 2 2**iT* , v 

El * [l =a a j -a -a -1 -a -a -j a a (4.14) 

P SI 


Crossmultiplication of the last twelve equations of (4.12), 
their addition and simplification of the expression gives the 
following result 
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-Vl V N -Yo V 2 V N -^oVn "Wn -Wj 

~ y o V 7 V N -y 0 V 8 V N ~ y o V 9 V N - y o V 10 V N - y o V ll V N - y o V 12 V * 

.7/-. , , , _* * *■ 2& » * * 

+12y o V N V N -y l V l E a “ y l a v 2 E a y l V 3 E a +jy l V 4 E a 
+ct2y l v 5 E a +0Iy l v 6 E t +y l v 7 E t + “* v 8 E a + “ 2 * y l v 9 E a 


- ^ y l V 10 E a - a2y l v ll E a _ccy ! V 12 E a + 12y l E a E a 


" S a + S b + ••• + Si 


(4.15) 


Noting that the first thirteen terms of eqn. (4.15) are V* 
times 1^ and the injected neutral current 1^ = 0, we can 
write eqn. (4.15) as 


£tc y O^* o - 

“ y l ~ a y l “ a y l jy l a y l ay l y l a y l a y l “JYi “ a Yl ~ a Yi j 


* * * 
S a +S b + ... + s x 


(4.16) 


If the neutral is earthed through an admittance y^ Q , eqn. (4.4) 
is modified to 


*N " I 13 +I 14 + * * * + *24 + y No V N 


(4.17) 
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Thus, from eqn. (4.13), we obtain 


I N = t-Y 0 “ y o ~ y o “ y o ” y o “ y o ** y o ~ Y o " y o ~ y o ~ Y o ~ Y o 12y o + 

(4.18) 

■if* 

Substituting the value of from eqn. (4.18) in eqn. (4,15), 

we get with l N = 0, the following 


+ y Nol 


V, 


12 


V. 


N 


r * 2* . 2 

L-Yi -<* Yi Y X 3Y 1 a Y 1 a Yi Y x 


*- 2 * 2 
a Y X « Yi -JYx “« Yi 





y No^ V N 


(4.19) 


where £S = S + + 

a b 


* 


+ S n 


Neglecting the second term in the R.H.S. of eqn. (4.19), the 
final equation of a twelve phase synchronous generator with all 
the types of unbalances included in it is obtained by 
collecting the first twelve equations of (4.12), eqn. (4.18) 
and eqn. (4,19). The concise form of the final eqn. is 


m ** 

v 12 

Yp 

-h - 7 o 


12 ~~ 
V 


12 

X P 

Y 2 

12y^ 0 


E a 

= 

* . * 
ES /B, 

^.yT 

1 0 

0 12 V Y No_ 


_ V N _ 


0 

'"BN** 


(4.20) 
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where, 

Y x = =[l a a 2 j -a 2 * -a* -1 -a -a 2 -j a 2 * a*] T y x , 

Y 2 = [1 a' a 2 * -j -a 2 -a -1 -a* -a 2 * j a 2 a] T y x , 

Y = [11111111111 l] T y 

o 

where a = e j2n/l2 _ 1 ^ 30 ° . 

Eqn. (4.20) is the complete eqn. of a twelve phase machine 
capable of including all the unbalances present in it. 

4.3 MODELLING OF TWELVE PHASE TRANSFORMERS 

Transformers used for multiphase systems are required for 
obtaining higher phase order conversion from a given three 
phase system. A twelve phase conversion can be obtained by 
using wye/star or delta/star three phase/twelve phase trans- 
formers. The feasibility of a twelve phase transformer and 
the procedure to obtain twelve phase connection from a three 
phase transformer has been given by Uma Pal et al. [24]. 

There may be other ways of obtaining twelve phase conversion 
also but in the present work, the models of the two types 
of transformers mentioned above have been developed. 
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The transformers are rep res 

in series with the equivalent p ed as ideal transformers 
windings. Consider a wye/star * le l ka 9 e impedances of the 
phase transformer shown in Fin ected » three phase/twelve 
former is represented by an itl * * A nominal r ^tio trans- 
the p.u. leakage impedance 2 of ^ ansforraer in series with 
and current relations at the + Q 6 Wlndin 9 s * The voltage 
mer are : s 0:f the ideal transfor- 


v “ = ; 
!p? = [m]i 3 ; 


where. 


V 3 ^ 1 r 

P " f [L] v 1 ? 

P 

i 3 i 

P f tL] & 

P f 


V P 

v “ - [Vi v 2 ... V 12 ] T , 

*p = f 1 . *b 

xj 2 = [Xj. I 2 ... i 12 ] t , 


(12x12) order diagonal 


and 


1 1 0 0 0 -1 , 


0 0 0 1 1 1 

0 - 1 - 1-10 0 


ra *tri x = [ y ]-l 

1 *’ 1 0 0 0 1 

0 0 -1 -1 -1 

0 i . 

1 1 0 0 


0 


( 4 . 21 ) 

( 4 . 22 ) 
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The terminal voltages of the transformer in terms of the 

12 

induced voltage vector V on the 12-phase side is given by 

V p 2 = v£ 2 - [Z] I^ 2 (4.23) 


The matrix L is given by 



2 

1 -1 

0 

1 

-1 -2 - 

-1 

1 

0 

-1 1 

L = 

1 

0 -1 

1 

2 

1 -1 

0 

1 

-1 

-2 -1 


1 

-1 -2 ■ 

-1 

1 

0 

1 

1 

2 

1 

-1 0 

From eqns. 

(4.21) 

and 

(4.23), we 

obtain 



Ip 2 = [y] [M]VJ* - [y]v£ 2 (4.24) 

Also from eqns. (4.22) and (4.24), we get 

Ip = £ [L] y[M]V^ - | [L] y V* 2 (4.25) 


Thus the nodal representation of a wye/star, three phase/ 
twelve phase transformer may be written as follows : 



[L] [y][M] 
-[y][M] 


- £ [L] [y] 
[y] 



(4.26) 



The (15x15) order nodal admittance matrix of such a transformer 
may thus be obtained by expanding eqm (4.26). 



In a similar manner, we can obtain the nodal admittance 
natrix of a delta/star, three phase/ twelve phase transformer 
shown in Fig. 4,3. Considering the turns ratio of such a 
transformer to be a-j. : 1, its general equation may be written 
as 


I 3 

p 


k 1 [L3[y3 [J] 

-k 2 [L][y] 


V 3 

p 

T 12 

* 

-k 3 [y]CJ] 

[y] 


v 12 
_ p 


The value of the constants are : 


For a wye/star transformer : 

k l “ JJT 5 *2 = 8o 7“ * k 3 “ * J = M 

For a delta/star transformer : 





k = -I— ; J = K 
f3a T 


where, 


K J 


ill 1000 -1 -1 -1 000 

-1 0 0 0 1 1 1 0 0 0 -1 -I 

0 0 -1 -1 -1 0 0 0 1 1 1 0 


Thus we see that the transformer models obtained from eqn. 
(4.27) are quite simple and reasonably accurate. In a similar 

obtain the models for other types of transformer 


manner, we can 




SCHEMATfC REPRESENTATION OF A THREE PHASE /TWELVE PHASE 
DELTA/STAR TRANSFORMER WITH TURNS RATIO o^T * 1 




4.4 MODELLING OF TWELVE PHASE TRANSMISSION LINES 


The representation of twelve phase transmission lines in 
terms of phase impedance matrix, equivalent 7t-circuit and ABCD 
parameters has been developed. Three phase equivalent represen- 
tation of a twelve phase line and vice versa has also been given. 

a) Phase Impedance Matrix Representation 

A twelve phase short untransposed transmission line posse- 
ssing cyclic symmetry can be described by its phase impedance 
matrix as below ; 


7 12 
^P 



(4.28) 


The matrix can be fully diagonalized by using symmetrical 
component twelve phase transformation. However, eqn. (4.28) 
cannot represent medium and long lines. 



b) it-Circuit -Representation 

The three phase and six phase nominal equivalent it-circuit 
representation can be extended to twelve phase systems as well 
in order to represent medium and long lines as shown s 


XV*J 




Eqrw (4.29) can be derived by employing the symmetrical lattice 
equivalent circuit of the twelve phase transmission line as was 
done for three phase lines. In eqn. (4.29) , the second sub- 
script S and R denote the sending end and receiving end res- 
pectively. 

c) ABCD Parameters Representation 

The twelve phase lines based upon ABCD parameters can be 
represented as below 



where A ,B ,C and D are square matrices of order 12. 

4.5 EQUIVALENT THREE PHASE REPRESENTATION OF A TWELVE PHASE 
TRANSMISSION LINE 


To carry out the analysis of a composite throe phase and j 

I 

twelve phase system, equivalent three phase representation of j 
the twelve phase part of it is required. For example, consideij 
a twelve-phase transmission line connected between throe-phase i 

I 

buses S and R via throe phase/twelve phase transformers as f 
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shovrn in Fig. 4.4. The phasor equation of the twelve-phase 
transmission line connected between buses S' and R’ in the 
impedance form is, 


,,12 ,JL2 _ 7 12 t 12 

v ^» V R* ^p x p 


S’ R 
,„12 


(4.31) 


12 

where Zp is the (12x12) phase impedance matrix of the 
twelve-phase transmission line. 


Let the transformers and T 2 be wye/star, three-phase/ 
twelve-phase of nominal turns ratio. Then, using the voltage 
relationships of eqn. (4.21), the three phase voltages $t the 

bus S and R may be expressed in terms of the twelve phase 

12 12 12 
voltages Vg, and Vp^, and the current Ip as shown. 


v i 

1 r ,,12 

- 3 L v s« 

+ k z t*i] if 

(4.32) 

< 

= 1 L V 12 
8 L V R* 

- | L [z 2 ] I* 2 

(4.33) 


where and [z 2 ] are the (12x12) diagonal matrices repre- 
senting the leakage impedances of the transformers T-^ and T 2 
respectively. 


Subtracting eqn. (4.33) from eqn. (4.32), we obtain 

V I‘ V R = 5 l C v I'- v R^ + 5 L if 

= l LZ p 2l p 2+ I L [[* 1 >t2 2 ]j if 


(4.34) 



Now using eqn* (4.22), we obtain 


v ! - ^ = I L + z i + z 2 ] M J p 

Therefore i Z^ eq = § L[z£ 2 + + z 2 ]M 

The three-phase equivalent admittance matrix of a 
line can either be obtained by inverting or 

get 

Vp >eq = [U+ 5 L yJ 2 [[z^+CZj]] m]- 1 [§ LYf m] 

where U is a (3x3) identity matrix. 

In a similar, we can obtain the three phase equivalents 
when other types of transformers are used. 

4.6 EQUIVALENT TOELVE-PHASE REPRESENTATION OF A THREE-PHASE 
LINE 

In a composite three-phase and twelve-phase system when 
the interest of investigation lies in the twelve-phase part of 
the network, it may be desirable to analyse the system entirely 
on twelve phase basis. For this purpose, a twelve-phase equi- 
valent of three-phase line is derived. For example, consider a 
three-phase line connected to the twelve phase buses via 
twelve-phase/three phase transformers as shown in Fig. 4.5. Tt 
three-phase lino is represented by 

v 3 ^3 ^3 j3 
V S' V R ! HP P 


(4.35) 

(4.36) 

twelve-phase 
derived to 

(4.37) 


(4.38) 





Using relationships (4.21) - (4.23), the voltage drop 
12 12 

Vg - Vj^ is expressed as 

V S 2 - V R 2 = Cs M 4 h + l Z l +2 4 ^P 2 (4 ' 39) 

Thus, from eqn. (4.39), the twelve phase equivalent representa- 
tion of the three phase line is given by 

Z J 2 eq = I M 4 1 + [Z 1 +Z 2 ] (4,40) 

Similarly, the twelve-phase equivalent admittance matrix is 

12 

either obtained by inverting Zp or derived in a similar 
manner to get 

Y p 2 oq =["*5“ V p l £l +z 2 j l' 1 £5 M Y p W.41) 

4.7 MODELLING OF TWELVE PHASE LOADS 

Twelve-phase loads may be modelled on the basis of constant 
impedance/admittance to ground. In this case also, we will j 

determine the equivalent three phase representation of a twelve- l 

| 

phase load and vice versa depending upon whether our interest 1 

I 

of investigation lies in the three-phase part or the twelve phase; 

I 

part of a composite three-phase, twelve phase system. [ 

- 1 

Consider a twelve phase load connected via a three-phase/ j 

I 

twelve-phase, wye/star transformer to the throo-phase bus as j 

j 

shown in Fig. 4.6. For the twelve-phase load, we can write j 


I 


XU? 


r 12 _ _12 t 12 

P " Z L X p 


We can 


or 


or 


3 

also write the volt, age Vp as 
^ - K + | L [ Z ] if 

Vp = 5 LZ L 2 MIp + | L[z] MI 2 
Z L,eq “ I L + Z] M 


(4.42) 


(4.43) 


(4.44) 


(4.45) 


Similarly, we can derive the equivalent three-phase load 
admittance as 


Y? = [U+ g h yJ 2 [zJm]" 1 [g L yJ 2 M] (4.46) 

Now consider a three-phase load connected to a twelve-phase 
bus through a twelve-phase/three-phase, star/wye transformer 
shown in Fig. 4.7. The load can be represented by constant 
impedance to ground as 

= Z? I 3 (4.47) 

p L p 

The voltage VT can be expressed as 

v p 2 = 5 M A. LI p 2 + ^ I p 2 


(4.48) 
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FIG. 4.4 A TWELVE PHASE TRANSMISSION LINE CONNECTED 

BETWEEN TWO THREE PHASE BUSES VIA THREE PHASE/ 
TWELVE PHASE TRANSFORMERS Tj AND T 2 
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FIG. 4.5 A THREE PHASE TRANSMISSION LINE CONNECTED 

BETWEEN TWO TWELVE PHASE BUSES THROUGH TWELVE 
PHASE/ THREE PHASE TRANSFORMERS T, AND T 2 


Ill 



Three Phase 3-0/12-0 Twelve Phase 

Bus Transformer Load 

FIG. 4.6 A TWELE PHASE LOAD CONNECTED TO THREE 
PHASE BUSBAR. 



Twelve Phase 12-0/3-0 Three Phase 

Bus Transformer Load 


FIG. 4.7 A THREE PHASE LOAD CONNECTED TO TWELVE 
PHASE BUSBAR 
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Hence the twelve phase equivalent impedance is 

z“ eq = £ M Zj* L + Z (4.49) 

Similarly, the twelve-phase equivalent admittance may be 
derived as 

= I U + 5 M Y L LZl" 1 [$ M Yl L] (4.50) 

4.8 CONCLUSIONS 

Although multiphase systems are only of theoretical 
interest today, twelve-phase transmission system may be 
feasible if double circuit six phase lines come into existence, j 

This is quite likely to happen in future because of the various J 

' )| 

advantages of six-phase systems over three-phase systems. In \ 
this chapter, the mathematical modelling of various twelve-phas<j 
elements such as twelve phase machines, transformers and trans- 
mission lines has been given. This modelling may be employed t 
study the feasibility of twelve-phase systems. 

I 

t 



CHAPTER 5 


CONCLUSIONS 

In the present work, an attempt has been made to study 
some aspects of mathematical modelling and analysis of multi- 
phase power systems. Although, multiphase systems are of 

r ** 

theoretical interest today, studies have revealed that they 
may become a reality in future. To get an insight of the 
behaviour of the multiphase systems, it is necessary to develop 
the mathematical models of the different components present in 
it so that detailed analysis like load flow analysis, fault 
analysis, transient stability analysis, etc. may be conducted. 

j 

It has been stressed, that, the mathematical modelling j 

of the elements of a power system in the phase frame of I 

reference gives a realistic picture of the system as it takes j 

f 

into account the inherent unbalances present in it. The j 

detailed models of three phase machines, transformers and 

transmission lines based upon their symmetrical lattice equi- , 

f 

valent circuits have been given. j 

As a direct extension of three phase mathematical mode- : 
lling, the representation of elements of multiphase (six : 

I 

phase and twelve phase) power systems in phasor coordinates f 
are given. The equation of the multiphase generators to ; 

include the effects of unbalances has been derived in the j 
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Matrix form. The multiphase transformers have been modelled 
based upon their symmetrical lattice equivalent circuits and 
their nodal relationships. The representation of the trans- 
mission line based upon its phase impedance matrix, -rt-equiva- 
lent circuit and ABC -parameters have been given. The multi- 
phase loads are modelled as constant admittances to ground. 
Equivalent three phase representations of the multiphase 
transmission lines and vice versa have also been developed. 

To show the advantages of the modelling in phasor 
coordinates, the load flow analysis of a three phase power 

system network and a mixed three-phase, six-phase power I 

! 

system network was carried out. With unbalances introduced ( 

i 

t 

in the three phase power system, it was found that there is j 
a considerable unbalance in the voltage magnitudes and power j 
flows of different phases. The effect of converting an | 

existing double circuits three phase line to a six phase j 

line was studied. It was found, that, the introduction of j 

six phase lines in the system improves the voltage magnitudes | 

! 

at the buses, reduces its angle increases the real power flow 
and reduces the reactive power flows from the generators. It 
was also observed, that, the six phase lines tend to reduce ’ 
slightly the degree of unbalance in the voltage magnitudes of j 
different phases at a bus. | 


I 



Although sufficient work seems to have been done on the 
theoretical feasibility of six phase transmission systems, 
more realistic picture of six phase systems can be had by 
conducting practical tests on installed six phase transmission 
lines and the results compared with the theoretical results. 
Further, work in this area may be pursued by taking a complete 
six phase system consisting of six phase generators, transfor- 
mers and transmission lines. The mathematical modelling 
developed for the twelve phase systems in the present work may 
be employed to carry out the detailed load flow studies, fault 
analysis and transient stability analysis of such systems. This 
will help in getting an overview of the practical feasibility 
of twelve phase transmission systems which may be a direct 
consequence of a double circuit six phase line. 
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